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Project  Aims: 

1.  To  develop  a  2-D  High  Performance  Computing  grid  point  model  of  the  ionosphere  in  the  region  of  the 
Tasman  International  Geospace  Environment  Radar  (TIGER)  coverage  area.  The  grid  point  model  will 
be  2-D  in  ground  range  and  height,  and  each  grid  point  will  specify  ionospheric  layer  parameters 
determined  from  real  TIGER  ground  scatter  data  using  a  new  inversion  technique  we  have  developed. 
This  modelling  effort  will  help  to  improve  our  understanding  of  ionospheric  structure  on  global  scales, 
including  gradients  in  electron  density,  the  occurrence  of  high-latitude  ionospheric  troughs  and  sporadic 
E,  and  the  associated  generation  of  ionospheric  irregularities. 

2.  To  extend  the  2-D  grid  point  model  to  a  3-D  grid  point  model  of  the  ionosphere  within  the  footprint  of 
TIGER.  TIGER  routinely  records  echoes  on  16  different  beam  directions  so  by  applying  the  2-D  grid 
point  approach  to  each  of  the  beams  a  3-D  description  of  the  ionosphere  may  be  achieved.  We  have  also 
developed  a  3-D  analytic  ray  tracing  technique  (not  yet  published)  that  can  be  applied  to  the  3-D  grid 
point  ionospheric  model  to  further  test  the  integrity  of  the  inversion  results. 


Performance  against  Projected  Milestones: 


First 

Milestones: 

Outcome 

Quarter 

•  Synthesize  TIGER  Backscatter  results  using  inversion  technique 

COMPLETED 

•  Develop  a  2-D  HPC  grid  point  model  for  the  E  layer 

COMPLETED 

Second 

Milestones: 

Quarter 

•  Develop  a  2-D  HPC  grid  point  model  for  the  FI  layer 

COMPLETED 

Third 

Milestones: 

Quarter 

•  Develop  a  2-D  HPC  grid  point  model  for  the  F2  layer 

COMPLETED 

•  Determine  the  accuracy  of  the  2-D  HPC  grid  point  model 

COMPLETED 

Fourth 

Milestones: 

Quarter 

•  Develop  a  3-D  HPC  grid  point  model 

COMPLETED  - 

PROTOTYPE 

DEVELOPED 

•  Synthesize  TIGER  Backscatter  results  using  3-D  HPC  SMART 

COMPLETED  - 

PROTOTYPE 

DEVELOPED 

•  Compare  the  synthesized  results  to  the  real  results 

COMPLETED - 

PROTOTYPE 

DEVELOPED 
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Outcomes: 

As  indicated  in  the  table  showing  performance  against  projected  milestones,  the  project  aims  have  been 
achieved.  In  particular  the  extension  of  the  analysis  techniques  to  three  dimensions  has  been  developed  and 
tested  in  prototype  form,  but  importantly  all  the  software  tools  that  we  set  out  to  develop  have  been 
developed  so  that  the  backscatter  inversion  techniques  can  now  be  applied  in  further  developments  such  as 
those  outlined  in  the  attached  white  page  on  a  proposed  follow-on  project. 


The  following  were  proposed  as  outcomes  in  the  work  Plan: 


Outcome 

Detail 

Conference  Presentations  and 
Papers 

Present  papers  at  one  national  and  two  international  conferences. 

Submit  at  least  one  paper  for  refereed  publication. 

Grants  Submitted  and/or 
Received 

Submit  at  least  one  application  for  government  grant  and  industry 
funding  for  TIGER  related  research 

The  following  were  the  actual  achievements: 

1.  Papers  (copies  attached  as  Appendices) 

1.1.  Dyson,  P.  L.,  M.  L.  Parkinson,  M.  Condeand  J.  C.  Devlin. 

New  Investigations  of  Magnetosphere-Ionosphere-Thermosphere  processes  and  their  Impact  on 
Space  Weather. 

5th  NSSA  National  Space  Science  Conference,  Melbourne,  Australia,  pp  1-17,  14-16  September, 

2005.  ISBN:  0  86459  3740.  (Refereed  Conference  Paper  National/International  Conference) 

1.2.  Norman  R.  J.  and  P.L.  Dyson  and  M.  L.  Parkinson. 

Inversion  of  HF  radar  backscatter  data  in  order  to  map  and  model  the  ionosphere. 

SuperDARN  2005  Workshop,  Ambleside,  UK,  16-20  May,  2005.  (International  Conference). 

1.3.  Norman  R.  J.  and  P.L.  Dyson  and  M.  L.  Parkinson. 

Mapping  the  Ionosphere  using  a  HF  Radar  Backscatter  Inversion  Technique. 

WARS06  (Workshop  on  Applications  in  Radio  Science),  Leura,  Australia,  8  pages,  15-17  February, 

2006.  (Accepted  for  publication  in  proceedings  as  a  refereed  paper.  WARS06  is  a  national 
conference). 

1.4.  Norman,  R.  J.  and  P.L.  Dyson. 

HF  Radar  backscatter  inversion  technique. 

Radio  Sci .,  (In  Press)  January  2006.  (Refereed  journal  publication). 


2. 


Grants  submitted 

2.1.  Title: 

Pis: 


Grant  Body: 
Status: 


A  consolidated  approach  to  ionospheric  studies  with  digital  radar  systems 

Dyson  and  Devlin,  La  Trobe  University,  Australia 

Lester,  Leicester  University  UK 

Bristow,  University  of  Alaska  Fairbanks,  USA 

Australian  Research  Council 

Unsuccessful 


2.2.  Title: 

Cl 

Grant  Body: 
Status: 


Investigations  of  Space  Weather  and  the  Mesosphere  using  the  TIGER  Radar 
Dyson,  La  Trobe  University,  Australia 
Australian  Antarctic  Science  Program 

Awarded  A$  17,000  for  2006/07.  Project  approved  for  5  years,  funding  to  be 
requested  on  an  annual  basis. 


In  summary,  the  proposed  outcomes  have  been  achieved. 


Prof  Peter  L.  Dyson 
Principle  Investigator 
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Appendix  1 


Dyson,  P.  L.,  M.  L.  Parkinson,  M.  Condeand  J.  C.  Devlin. 

New  Investigations  of  Magnetosphere-Ionosphere-Thermosphere  processes  and  their  Impact  on 
Space  Weather. 

5th  NSSA  National  Space  Science  Conference,  Melbourne,  Australia,  pp  1-17,  14-16  September, 
2005.  ISBN:  0  86459  3740.  (Refereed  Conference  Paper  National/International  Conference) 
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New  Investigations  of  Magnetosphere-Ionosphere- 
Thermosphere  processes  and  their  Impact  on  Space 

Weather 


12  2  2  3 

Peter  L.  Dyson  ’  ,  Murray  L.  Parkinson  ,  M.  Conde  and  J.  C.  Devlin 


1 CRC  for  Satellite  Systems 
2 Department  of  Physics 
department  of  Electronic  Engineering 
La  Trobe  University ,  VIC  3086 
AUSTRALIA 


Summary:  Space  weather  impacts  on  some  important  modem  technologies,  particularly 
those  that  depend  on  satellites. 

The  La  Trobe  Space  Physics  Research  Group  studies  the  basic  physical  causes  of  space 
weather  phenomena  in  geospace,  the  interface  between  Earth  and  interplanetary  space.  The 
investigations  are  based  primarily  on  three  experimental  techniques:  the  Tasman  International 
Geospace  Environment  Radar  (TIGER)  dual  HF  radar  system;  optical  spectrometers  used  to 
study  auroral  and  airglow  emissions;  and  GPS  tomography  using  the  Australian  satellite 
FedSat.  Future  plans  include  development  of  rocket  instmments  to  be  flown  from  Woomera. 

This  paper  will  first  present  a  brief  overview  of  some  relevant  aspects  of  solar  wind- 
magnetosphere-ionosphere  coupling  that  are  the  main  drivers  of  space  weather  in  geospace. 
The  La  Trobe  research  program  will  be  described  and  recent  research  results  of  storm-related 
stmcture  in  the  plasmasphere,  sub-auroral  motions  and  thermospheric  wind  fields  will  be 
presented  and  discussed.  Finally,  future  research  directions  will  be  described. 

Keywords:  space  weather,  aurora,  magnetic  storm,  geospace,  HF  radar,  SuperDARN, 
ionospheric  tomography,  Fabry-Perot  interferometer. 


Introduction 

Space  weather  refers  to  the  influences  some  important  conditions  on  the  Sun  and  in  the  solar 
wind  have  on  the  magnetosphere,  ionosphere  and  thermosphere.  This  solar-terrestrial 
environment,  also  called  geospace,  is  dominated  near  the  Earth  by  the  interaction  of  the 
particles,  fields,  and  radiation  environment  from  the  Sun  with  Earth’s  space  plasma 
environment  and  upper  atmosphere.  Emissions  from  the  Sun  impact  on  the  Earth’s  outer 
atmosphere  producing  a  variety  of  effects  including  large  scale  disturbances  in  the  Earth’s 
magnetic  field  and  the  aurora.  While  x-rays  and  high  energy  particles  generated  by  solar 
flares  travel  to  the  Earth  in  minutes,  other  emissions  such  as  coronal  mass  ejections  and 
features  in  the  solar  wind  take  up  to  several  days. 

Space  weather  impacts  on  some  important  modern  technologies  and  particularly  those  that 
depend  on  satellites:  long-distance  communications  (e.g.  broadcasts  TV,  radio,  telephone  and 
internet);  GPS  navigation  and  position  location;  remote  sensing  and  weather  forecasting. 
Energetic  particles  associated  with  space  weather  events  can  damage  satellite  electronics 


causing  satellite  malfunction  or  failure,  and  large  magnetic  field  fluctuations  and  related 
currents,  electric  fields  and  plasma  motions  in  the  ionosphere  and  magnetosphere  also 
produce  significant  effects,  such  as  disruption  of  radio  communications  and  over-the-horizon 
radar,  and  reduced  accuracy  in  GPS  navigation,  and  nanotechnology  lithography.  Many 
detailed  aspects  of  space  weather  phenomena  are  poorly  understood,  limiting  our  ability  to 
predict  their  occurrence  and  impact. 

The  solar  wind  plasma  carries  with  it  the  Interplanetary  Magnetic  Field  (IMF)  that  interacts 
strongly  with  Earth’s  magnetic  field,  confining  it  in  space,  particularly  on  the  dayside  where 
the  magnetopause  is  formed  at  the  pressure  balance  between  the  solar  wind  and  the  Earth’s 
magnetic  field  (Fig.  1).  Under  quiescent  solar  wind  conditions  the  magnetopause  forms  at  -10 
Re  (where  Re  is  the  radius  of  the  Earth,  -6300  km).  During  disturbed  conditions  the 
magnetopause  moves  closer  to  Earth,  to  -6Re,  or  even  closer  during  very  large  events. 

The  solar  wind  is  rarely  steady  for  long  periods  of  time  and  the  IMF  direction  can  change 
back  and  forth  between  having  a  northward  component  and  southward  component,  as  shown 
in  Fig.  1.  When  the  IMF  is  northward,  its  direction  perpendicular  to  the  equatorial  plane  on 
the  dayside  of  the  Earth  is  the  same  as  the  direction  of  the  terrestrial  magnetic  field.  However, 
when  the  IMF  is  southward  the  two  field  directions  are  opposite  and  IMF  field  lines  can 
connect  directly  to  terrestrial  field  lines,  dramatically  increasing  the  coupling  of  energy  from 
the  solar  wind  into  the  magnetosphere  and  consequently  into  the  lower  ionosphere  and 
thermosphere.  The  coupling  is  not  restricted  to  magnetic  energy;  it  includes  the  coupling  of 
solar  wind  plasma  into  the  terrestrial  environment.  Of  course  the  interaction  is  not  restricted 
to  the  dayside  as  is  shown  in  Fig.  1,  which  only  shows  a  cross-section  perpendicular  to  the 
ecliptic.  The  solar  wind-magnetosphere  interaction  occurs  over  the  three  spatial  dimensions 
and  Fig.  2  shows  that  the  result  is  a  complex  system  characterized  by  different  regions  of 
fields,  currents  and  plasma. 


tnferptanetary  Fi'e/d  Southward 


Interplanetary  Fietd  Northward 


Fig.  1:  The  interaction  between  the  Interplanetary  Magnetic  Field  (IMF)  and  Earth ’s  mag¬ 
netic  field  in  a  plane  perpendicular  to  the  ecliptic  (after  [1]). 
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Fig.  2:  Schematic  of  the  magnetosphere  showing  currents,  fields  and  plasma  regions  (after 

m 

Of  specific  interest  here  is  that  the  field  lines  throughout  the  vast  region  of  the  magnetosphere 
that  is  most  affected,  connect  back  to  high  latitudes  on  the  Earth’s  surface.  As  a  result,  at 
heights  of  90  -  1000  km,  it  is  the  high  latitude  ionosphere  and  thermosphere  that  are  affected 
most  directly  by  the  solar  wind’s  interaction  with  Earth,  and  therefore  it  is  at  these  latitudes 
that  many  of  the  near-earth  space  weather  systems,  such  as  the  aurora,  generally  occur.  By 
studying  these  high  latitude  regions  we  learn  not  just  about  their  response  to  changing  space 
weather,  but  about  the  processes  that  transfer  energy  from  the  solar  wind  to  the  Earth’s 
environment.  This  is  essential  information  for  space  weather  prediction  services. 
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Fig.  3:  Representation  of  ionospheric  electric  fields  in  the  northern  hemisphere  polar  cap  and 
auroral  zone,  as  well  as  the  resultant  plasma  flow  (after  [2]). 


A  major  feature  of  the  high  latitude  ionosphere  is  the  development  of  convection  cells  in  the 
ionosphere  in  which  the  ionospheric  plasma  is  transported  across  the  poles  from  the  dayside 
to  the  nightside  with  the  return  flow  at  lower  latitudes  (Fig.  3).  The  convection  comes  about 
through  E  x  B  drift  of  the  plasma  caused  by  a  dawn-to-dusk  electric  field  set  up  across  the 
magnetosphere  by  the  solar  wind  interaction.  This  field  is  conducted  down  the  magnetic  field 
lines  to  the  ionosphere,  and  coupled  with  the  almost  vertical  magnetic  field  at  high  latitudes, 
produces  significant  horizontal  E  x  B  convection  in  which  the  plasma,  which  in  the 
ionospheric  F-region  at  -300  km  altitude,  reaches  speeds  as  high  as  1-2  km/s.  Changes  in  the 
IMF  can  produce  large  scale  changes  in  ionospheric  convection  that  occur  within  minutes  of 
the  IMF  changes  reaching  the  magnetopause.  Collisions  between  the  ions  and  neutrals 
produce  similar  flow  structures  in  the  neutral  thermosphere  above  100  km  which  are 
superimposed  on  the  tidal  flows  produced  by  solar  heating  of  the  atmosphere. 

Many  other  space  weather  processes  occur  in  the  thermosphere  and  ionosphere,  including 
particle  precipitation  that  produces  the  aurora,  electrojet  currents  that  flow  in  the  E-region  of 
the  ionosphere  and  produce  magnetic  disturbances  on  the  ground.  These  disturbances  can  also 
trigger  atmospheric  gravity  waves  that  propagate  to  low  latitudes;  field-aligned  currents  that 
flow  between  the  magnetosphere  and  the  ionosphere;  and  flow  channels  that  form  as  part  of 
the  detailed  structure  of  the  convection  flows. 

The  aurora  is  visually  spectacular  but  can  also  be  observed  at  wavelengths  outside  the  visible 
region.  Fig.  4  is  an  observation  in  the  far  ultraviolet  made  by  the  IMAGE  satellite.  The  aurora 
is  seen  to  form  an  oval  around  the  magnetic  pole,  varying  in  intensity  and  latitudinal  extent 
with  local  time.  The  ionospheric  convection  pattern  overlays  the  oval  and  the  enclosed  polar 
region  and  much  current  research  is  aimed  determining  and  understanding  the  details  of  the 
correspondence  between  the  aurora,  the  convection  and  other  space  weather  processes,  such 
as  those  listed  in  the  previous  paragraph. 


Fig.  4  The  aurora  australis  photographed  by  the  Far  UV Imager  on  the  IMAGE  satellite 
(http://son.nasa.gov/tass/auroralfriends/sat  image. him) 


To  understand  these  processes  ultimately  requires  coordinated  studies  using  a  world- wide 
network  of  satellite  and  ground  based  techniques  that  invariably  involve  international 
collaboration.  Obviously  instruments  that  measure  details  of  the  flow  structures  in  the  plasma 
and  neutral  atmosphere  at  high  latitudes  will  play  key  roles  in  addressing  the  gaps  in  our 
understanding  of  space  weather  processes  because  they  can  be  used  to  investigate  the 
development  and  subsequent  decay  of  many  of  the  phenomena.  In  this  paper  we  briefly 
describe  the  instrumental  techniques  being  used  and  proposed  by  the  Space  Physics  Research 
Group  at  La  Trobe  University  in  our  studies  of  geospace.  Some  recent  results  are  also 
presented. 


The  Tasman  International  Geospace  Environment  Radar  (TIGER) 

TIGER  is  a  dual  HF  radar  system  of  the  Super  Dual  Auroral  Radar  Network  (SuperDARN) 
type  [3],  [4]  and  in  fact  is  a  component  of  the  SuperDARN  network  of  radars  [5]  that  partially 
ring  the  north  and  south  polar/high  latitude  ionospheres  (Fig.  5). 


Fig.  5  The  SuperDARN  network:  southern  hemisphere  (left);  northern  hemisphere  (right), 
(after  R.  Barnes  http://superdarn.jhuapl.edu/index.html) 


Fig.  6  TIGER  locations  and  radar  footprints  (Courtesy  R.  Healy). 

One  TIGER  radar  is  located  on  Bruny  Island,  Tasmania,  the  other  near  Invercargill,  New 
Zealand  (Fig.  6).  An  aspect  of  the  high  latitude  phenomena  of  interest  is  that  they  produce 
irregularities  or  turbulence  in  the  ionosphere  with  structures  that  are  aligned  along  the  Earth’s 
magnetic  field  lines  and  which  scatter  HF  radio  waves.  Strong  backscatter,  which  sends 
signals  back  to  the  radar  site,  occurs  when  the  HF  radio  signals  become  perpendicular  to  the 
magnetic  field  lines  (Fig.  7). 


Fig.  7  Ray  paths  showing  ionospheric  backscatter  from  irregularities  aligned  along  the 
Earth ’s  magnetic  field  and  from  the  Earth ’s  surface. 


Only  a  small  fraction  of  the  radio  wave  energy  is  backscattered  and  at  some  frequencies  and 
propagation  angles  the  forward  propagating  wave  is  reflected  by  the  ionosphere  and  returns  to 
the  ground  where  it  is  reflected  and  so  propagates  to  even  greater  ranges  via  a  second 
ionospheric  ‘hop’.  As  a  consequence  ionospheric  backscatter  can  actually  be  detected  from 
0.5  and  1 .5  hops  (and  in  principle  from  even  greater  hops  depending  on  the  strength  of  the 
backscatter  and  sensitivity  of  the  radar).  Of  course  some  energy  is  backscattered  at  each 
ground  reflection  and  so  backscattered  signals  from  the  ground,  or  ocean  in  the  case  of  the 
TIGER  radars,  is  also  observed.  Another  source  of  echoes  is  backscatter  from  meteor  echoes 
at  ~  100  km  altitude. 

The  transmitting  antenna  for  the  radars  consists  of  a  phased  array  made  up  of  16  log-periodic 
antennas  placed  on  the  top  of  15  m  high  towers  (Fig.  8).  These  provide  a  beam  that  is  -3.6°  in 
azimuth  but  broad  in  elevation  so  that  rays  launched  at  different  elevations  effectively 
floodlight  the  ionosphere  out  to  distances  of  -3000  km.  The  beam  is  switched  sequentially 
through  16  steps  in  azimuth  giving  an  azimuthal  coverage  of  -52°  as  shown  in  Fig.  6.  Echoes 
are  received  on  both  the  transmitting  array  antenna  and  a  second  small  array  of  4  similar  log- 
periodic  antennas  placed  behind  the  transmitting  array  so  that  the  two  arrays  form  an 
interferometer  used  to  determine  echo  elevation  angle  (Fig.  8). 


Fig.  8  TIGER  Unwin  Antennas.  16  element  transmit/receive  antenna  (backround)  and  4 
element  receive  only  antenna  (forground)  [< Courtesy  M.  . 

The  standard  mode  of  operation  is  for  each  of  the  two  radars  to  send  out  test  transmissions 
over  8-20  MHz  and  select  the  frequency  that  has  the  greatest  range  extent  of  ionospheric 
scatter  and  then  for  the  two  radars  to  sweep  in  azimuth  in  a  synchronised  manner,  providing 
complete  scans  every  1-2  minutes.  The  radar  operating  system  is  very  flexible  and  other 
modes  can  be  programmed  which,  say,  vary  the  transmitting  frequency,  the  order  or  number 
of  azimuthal  beams  scanned,  and/or  the  time  for  a  complete  azimuthal  scan. 


TIGER  provides  a  unique  contribution  to  SuperDARN.  It  is  the  first  SD  radar  pair  to  be  located 
at  comparatively  low,  sub-auroral  latitudes,  allowing  new  studies  of  these  regions  that  are  critically 
important  to  understanding  what  is  called  the  magnetic  sub-storm  phenomenon,  when  the  Earth’s 
magnetic  field  in  the  midnight  sector  is  disturbed  and  space  weather  features,  such  as  the  aurora, 
move  equatorward,  often  beyond  the  field  of  view  of  other  SuperDARN  radars.  The  generation  of 
related  disturbances,  such  as  atmospheric  gravity  waves,  and  how  they  propagate  equatorward  over 
Australia  are  important  for  regional  communications  and  Australia’s  surveillance  operations  using 
the  Jindalee  Operational  Radar  Network  (JORN). 

Furthermore,  the  second  TIGER  radar  at  Invercargill,  named  Unwin  after  a  New  Zealand  scientist 
who  was  a  pioneer  in  the  development  of  auroral  radar,  is  one  of  the  few  SuperDARN  radars  that 
has  a  stereo  capability  enabling  it  to  operate  on  two  frequencies  simultaneously.  Thus  it  is 
effectively  two  radars,  and  one  channel  can  be  synchronised  with  the  Bruny  Island  (BI)  radar,  while 
the  other  channel  is  used  to  conduct  special  campaigns. 

For  each  echo  the  radars  record  the  range,  angle  of  arrival,  power,  line-of-sight  velocity 
(determined  from  the  echo  Doppler  shift)  and  Doppler  spectral  width,  usually  out  to  a  range 
of  3000  km.  Ionospheric  scatter  echoes  are  mostly  distinguished  from  ground  or  sea  scatter  by 
their  larger  velocities  and/or  spectral  widths.  Fig.  9  presents  a  BI  record  obtained  using  one 
beam  direction.  It  shows  the  spectral  width  of  echoes  as  a  function  of  radar  frequency.  It 
shows  two  similar  features  displaced  in  range  and  characterised  by  different  spectral  widths 
so  that  one  appears  black  and  the  other  red.  Each  feature  consists  firstly  of  a  broad  region  of 
echoes  for  which  the  range  hardly  changes  as  the  radar  frequency  is  changed.  These  are  due  to 
propagation  via  what  is  known  as  a  sporadic-E  layer,  at  about  100  km  altitude.  The  other 
characteristic  of  each  feature  is  a  broad  echo  trace  at  longer  ranges  that  increases  in  range  as 
the  radar  frequency  is  increased.  These  echoes  arise  from  propagation  via  the  F-region  that 
usually  has  its  peak  electron  density  at  -300  km  altitude.  The  low  spectral  width  echoes 
(black)  are  due  to  1-hop  sea  scatter  and  the  large  spectral  width  echoes  (red)  are  due  1.5-hop 
ionospheric  scatter. 

A  recent  discovery  made  with  TIGER  is  the  occurrence  of  Auroral  Westward  Flow  Channels 
(AWFC)  near  the  equatorward  edge  of  the  auroral  oval  (Parkinson  et  al.,  2003).  Fig.  1 1  shows 
observations  of  such  an  event  by  the  BI  radar.  In  the  radar  scan  shown  on  the  left  of  the 
figure,  the  echo  region  from  0.5  hop  propagation  shows  very  large  line  of  sight  velocities 
approaching  the  radar  on  the  eastern  side  of  the  scan  and  receding  to  the  west  on  the  western 
side.  These  line  of  sight  components  are  consistent  with  a  very  high  speed  left  flow  of 
ionisation  westward  across  the  radar’s  field  of  view.  The  orange  and  red  colours  on  the 
eastern  side  of  the  scan  show  that  there  is  latitudinal  structure  in  the  strength  of  the  flow, 
indicating  that  the  flow  consists  of  several  high  speed  channels  within  a  general  high  speed 
flow. 

The  plot  on  the  right  of  Fig.  1 1  shows  the  time  evolution  of  the  AWFC  as  it  appears  along  one 
radar  beam  direction.  The  flow  begins  very  abruptly  and  last  for  several  hours,  although  it 
reduces  in  speed  after  about  an  hour.  Magnetometer  data  from  Macquarie  Island,  which  is 
under  the  Unwin  footprint,  shows  that  the  onset  of  these  AWFC  channels  closely  coincides 
with  the  onset  of  magnetic  sub-storms  [7],  In  another  study  using  the  King  Salmon  SD  radar 
in  Alaska,  we  have  shown  that  AWFCs  form  in  both  hemispheres  (although  with 
asymmetries)  as  part  of  the  sub-storm  phenomenon  [8]. 


Fig.  9  Observation  of  1-hop  sea-scatter  and  1.5-hop  ionospheric  scatter.  Echo  range  as  a 
function  of  radar  frequency  colour  coded  according  to  Doppler  spectral  width. 
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Fig.  10.  TIGER  Bruny  Island  radar  scan  illustrating  the  use  of  echo  velocity  and  spectral 
width  to  identify  different  echoing  regions. 


Woe**  (m/*) 


Velocity  (m/s) 

-60  000  060  120  180  240  300  360  420  460 


10  11  12  13  14 


Fig.  11.  Observation  of  an  Auroral  Westward  Flow  Channel.  Radar  scan  (left)  showing 
strong  structured  westward  flow  in  0.5 -hop  ionospheric  scatter.  Single  beam  observations 
(right)  showing  sharp  onset  of AWFC  and  duration. 


The  existence  of  these  AWFCs  raises  a  number  of  important  questions  which  we  are  currently 
investigating: 

•  What  fraction  of  the  total  energy  dissipation  of  a  substorm  is  associated  with  AWFCs? 

•  How  closely  synchronised  are  the  start  of  AWFCs  and  magnetospheric  substorms? 

•  Can  AWFCs  precede  the  start  of  substorms,  and  act  as  a  trigger  for  them? 

•  Do  AWFCs  occur  under  geomagnetic  quiet  conditions? 

•  Can  we  resolve  the  multiple,  filamentary  convection  structure  embedded  within 
AWFCs,  and  how  large  are  the  peak  velocities  within  these  filaments? 

•  What  is  the  relationship  between  AWFCs,  radiation  belt  particles,  and  sub-visual 
auroral  red  (SAR)  arcs? 

•  What  plasma  instabilities  act  within  AWFCs,  and  are  AWFCs  a  source  of  atmospheric 
gravity  waves? 

Other  phenomena  studied  using  TIGER  include: 

•  Propagation  of  atmospheric  gravity  waves  from  the  auroral  region  [9] 

•  ULF  waves  and  micropulsations  [10],  [11] 

•  Signatures  of  magnetospheric  boundaries  in  HF  radar  data  [12],  [13] 

•  Mesospheric  winds  derived  from  meteor  echoes  [14] 

•  Sea- state. 

Future  areas  of  study  will  concentrate  on: 

•  Identifying  signatures  of  magnetosphere-ionosphere  coupling  and  magnetospheric 
boundaries. 

•  Auroral  westward  flow  channels  (AWFCs)  and  their  relationship  to  substorms. 

•  Ultra-low  frequency  (ULF ;  1-100  mHz)  wave  fields  imaged  by  the  TIGER  radars. 

•  The  propagation  of  ionospheric  disturbances  from  the  auroral  zone. 

•  Identifying  signatures  of  severe  space  weather  events. 

•  Developing  methods  of  determining  sea-state 

•  Identifying  tidal  motions  in  the  mesosphere-lower  thermosphere  from  the  analysis  of 
meteor  echoes. 

A  new  generation  digital  radar  is  being  designed  [15],  [16].  This  will  include  a  receiver  on 
each  antenna  as  opposed  to  the  current  two  receivers  -  one  for  each  antenna  array.  This  will 
provide  greater  flexibility  in  operation  permitting,  for  example,  beam  formation  in  post 
processing,  so  that  specific  phenomena  can  be  examined  and  tracked  in  more  detail  once  their 
occurrence  has  been  identified 


Optical  Observations  of  Thermospheric  Winds  and  Temperatures 

The  aurora  is  not  the  only  visible  emission  from  the  upper  atmosphere.  Atomic  oxygen  in  the 
thermosphere  and  OH  in  mesosphere  emit  sub-visual  airglow  that  can  be  detected  from  the 
ground  with  appropriate  instrumentation.  Two  important  atmospheric  parameters  are 
temperature  and  wind  and,  if  the  height  of  the  emitting  O  and  OH  is  known,  the  temperature 
can  be  determined  from  the  emission  line  width  and  the  line-of-sight  wind  component  from 
the  Doppler  shift  of  the  emission  line  peak.  The  same  atomic  oxygen  emissions  occur  in 
aurora  and  airglow  so  changes  in  thermospheric  wind  and  temperature  can  be  measured  by 
observing  a  single  emission  line.  Our  research  program  is  directed  at  observations  of  these 
emission  lines,  providing  measurements  at  atmospheric  heights  primarily  in  the  range  90-250 
km. 


A  widely  used  technique  for  observing  airglow  and  auroral  emissions  is  to  use  high- 
throughput  (-150-mm  plate  aperture)  Fabry-Perot  Interferometers  (FPIs)  and  observe  several 
directions  in  sequence,  e.g.  north,  south,  east,  west,  all  at  30°  elevation,  and  the  zenith.  We 
currently  operate  FPIs  in  Antarctica  at  Davis  and,  at  mid-latitudes  at  Beveridge,  near 
Melbourne.  An  instrument  operated  for  over  a  sunspot  cycle  at  Mawson,  Antarctica  and  a 
new  improved  replacement  instrument  is  under  construction.  While  many  important  results 
have  been  obtained  on  the  dynamics  and  thermodynamics  of  the  thermosphere  and 
mesosphere  [17] -[22],  the  spatial  resolution  of  the  instruments  used  has  been  limited  (-800 
km)  and  will  be  greatly  improved  in  the  new  instrument. 

The  spatial  resolution  of  the  existing  instruments  is  less  than  desirable  because  significant 
variations  with  much  smaller  scales  than  -800  km  occur  in  the  high  latitude  thermosphere. 
Ground  and  space  based  instruments  observe  vertical  wind  speeds  up  to  100  m/s  or  more, 
wind  shears  that  reverse  the  horizontal  flow  within  less  than  100  km,  and  composition 
perturbations  that  are  equivalent  to  more  than  100  km  of  vertical  transport.  Currently  these 
features  are  too  small  to  be  resolved  easily  by  observations  or  explained  by  models.  These 
small  scales  cannot  be  ignored.  Local  mixing  and  intense  heat  deposition  are  positive-definite 
phenomena;  there  are  no  corresponding  processes  to  negate  their  effects  through  local  “un¬ 
mixing”  or  through  equally  intense  local  cooling.  Their  cumulative  impact  modifies  the  large- 
scale  structure  of  the  upper  atmosphere  and  ionosphere  at  altitudes  above  100  km,  including 
heights  of  low  Earth  orbiting  satellites. 

To  investigate  these  problems,  Conde  [23]  developed  an  imaging  FPI  located  in  Alaska  with  a 
resolution  of  -100  kms.  One  result  is  shown  in  Fig.  12  which  shows  the  collapse  of  the 
thermospheric  horizontal  wind  field  during  an  auroral  break-up  phase  over  Alaska.  This 
collapse  occurred  in  a  few  minutes  and  was  accompanied  by  a  large  upward  vertical  wind.  As 
well  as  our  interest  in  understanding  the  processes  in  the  thermosphere  during  these  events, 
we  are  also  interested  in  what  role  such  gradients  play  in  modifying  the  ionosphere. 


Fig.  12.  Observations  of  thermospheric  horizontal  winds  over  Alaska  obtained  using  an 
imaging  Fabry-Perot  Interferometer  [after  24]. 


Conde  is  now  developing  a  new  state  of  the  art,  hyper- spectral  imaging  Fabry-Perot  optical 
interferometer  (HI-FPI)  to  be  deployed  at  Mawson,  Antarctica  (Fig.  13).  The  all-sky  lens  is 
used  to  obtain  a  150°  view  of  the  sky.  The  interference  filter  has  a  pass  band  of  a  few 
angstroms  and  is  used  to  isolate  the  emission  line  of  interest.  A  filter  wheel  provides  a  set  of 
filters  so  that  different  emission  lines,  coming  from  different  heights  in  the  mesosphere- 
thermosphere  region  can  be  selected.  The  etalon  produces  the  fringes  which  are  detected  by 
the  CCD  camera.  The  separation  of  the  etalon  plates  is  scanned  piezoelectrically  through  a  set 
range  so  that  each  camera  pixel  effectively  acts  as  an  individual  FPI  giving  a  fringe  scan 
associated  with  a  specific  point  in  the  sky. 

Based  on  the  operation  of  the  Alaskan  system,  on  a  single  night  of  active  aurora,  the  new 
system  will  acquire  35,000  or  more  spectra  of  very  high  signal-to-noise  ratio.  By  contrast,  the 
simpler  FPI  that  operated  at  Mawson  during  the  1980’s  and  1990’s  typically  gathered  only 
10,000  to  15,000  noisier  spectra  over  an  entire  year.  Thermospheric  temperatures  are  of  the 
order  of  1000K  and  uncertainties  are  expected  to  be  less  than  10K  in  the  zenith,  and  less  than 
20K  near  the  edge  of  the  field.  Corresponding  worst-case  line-of-sight  wind  uncertainties  are 
estimated  to  be  3  m/s  and  6  m/s  respectively.  Vertical  winds  in  particular  may  be  close  to  zero 
at  times  and  the  largest  winds  are  over  150  m/s. 

We  installed  an  improved  FPI  at  Davis  over  the  2003-04  summer.  Its  basic  optical  elements 
are  similar  to  that  shown  if  Fig.  13  except  that  it  has  a  periscope  rather  that  an  all-sky  lens  and 
the  etalon  has  a  fixed  spacing.  In  this  system  a  single  CCD  image  refers  the  (limited)  region 
of  the  sky  defined  by  the  pointing  direction  of  the  periscope.  The  instrument  is  capable  of 
quite  high  time  resolution,  able  to  obtain  a  measurement  is  tens  of  seconds  compared  to 
minutes  for  the  instrument  it  replaced.  It  is  also  computer  controlled  and  largely  operated  as  a 
remote  instrument  from  La  Trobe  via  the  ethernet. 


Fig.  13.  Schematic  of  the  new  hyper-spectral  imaging  Fabry-Perot  optical  interferometer 

(HI-FPI). 


The  new  instrument  will  be  used  to  study  the  following  phenomena  and  it  will  be  supported 
by  the  current  FPI  instrument  at  Davis. 


Small-scale  (100  km  and  less)  phenomena:  These  relatively  unexplored  spatial  scales  are 
where  the  new  instrument  will  excel.  Of  course  there  is  little  point  observing  such  scales 
without  supporting  data  to  show  what  is  driving  them.  We  have  an  FPI  at  Davis  and  both 
Mawson  and  Davis  are  favourably  located  beneath  the  field-of-view  of  the  Syowa  East 
SuperDARN  (SD)  radar,  which  observes  ionospheric  convection.  This  will  enable  us  to 
address  questions  such  as: 

1.  What  is  the  response  time  of  the  neutral  atmosphere  to  changes  in  ionospheric 
convection  (as  measured  by  SD),  and  how  does  this  vary  with  altitude? 

2.  At  what  rate  must  mechanical  energy  be  supplied  to  sustain  observed  horizontal  wind 
shears? 

3.  Do  small-scale  space  weather  disturbances  penetrate  downward  from  the  F-region  into 
the  E-region,  or  are  they  excited  simultaneously  throughout  the  thermospheric  column? 

4.  What  are  the  relative  contributions  of  local  heating  and  horizontal  divergence  to  driving 
large  vertical  wind  events?  How  does  this  vary  with  altitude? 

5.  With  the  higher  spatial  resolution  and  smaller  uncertainties  of  the  proposed  instrument, 
are  we  able  to  observe  auroral  excitation  of  thermospheric  gravity  waves?  If  so,  can  we 
explain  why  they  appear  to  propagate  preferentially  poleward,  at  least  in  the  middle  and 
upper  thermosphere?  (eg  [25],  [26]). 

Regional-scale  (100-1000  km)  phenomena:  Greet  et  al.  [21]  showed  that  regional  structure 
(hundreds  of  km)  in  the  neutral  wind  field  could  be  resolved  by  combining  observations  from 
Mawson  and  Davis  (-650  km  apart).  With  the  new  instrument,  we  can  expand  this  work 
considerably.  Topics  of  interest  include: 

1.  How  do  the  regional-scale  wind  fields  compare  with  SD  convection  patterns? 

2.  Can  we  observe  disturbances  propagating  in  local  time  -  for  example  in  response  to  a 
westward  travelling  surge? 

3.  Is  there  a  measurable  zonal  gradient  of  F-region  neutral  temperature?  If  so,  under  what 
circumstances? 


Satellite  Tomography  of  the  Plasmasphere 

Soon  after  satellites  were  first  launched  into  Earth  orbit  it  was  recognised  that  satellite  radio 
transmissions  could  be  used  to  determine  the  Total  Electron  Content  (TEC)  along  the  ray 
paths  from  satellites  to  ground  receiving  stations.  Over  the  years  this  has  become  an  important 
means  of  studying  the  ionosphere  and  the  higher  altitude  plasmasphere.  Several  methods  can 
be  used  based  on  measurements  of  different  radio  wave  properties  such  as  Faraday  rotation, 
Doppler  shift,  and  the  time  delay  of  signal  propagation  (which  is  possible  with  GPS  systems). 
More  recently  it  was  recognised  that  observations  from,  say,  several  GPS  satellites,  can 
provide  a  mesh  of  intersecting  ray  paths  suitable  for  tomography  [27]  and  this  technique  has 
been  used  by  our  group  to  study  the  response  of  the  ionosphere  in  the  Australian  region  to 
space  weather  impacts  [28],  [29], 

In  principle,  the  use  of  GPS  satellites  for  tomography  provides  the  capability  of  determining 
the  3-D  structure  of  the  ionosphere-plasmasphere  system  out  to  -22,600  km,  the  orbital  height 
of  the  satellites.  However  in  reality,  the  plasma  density  falls  off  rapidly  above  and  below  the 
F-region  ionospheric  peak,  typically  at  250-300  km  altitude,  so  that  GPS  tomography  using 
ground  stations  provides  excellent  detail  of  the  F-region  of  the  ionosphere,  but  not  of  the 
plasmasphere  above. 


The  Australian  satellite  FedSat  (www.crcss.csiro.au/fedsat/),  orbits  at  an  altitude  of  ~800  km, 
well  above  the  F-region  peak  and  so  provides  and  excellent  platform  for  GPS  tomography  of 
the  plasmasphere  which  is  a  region  about  which  little  is  known  of  its  detailed  structure  and 
response  to  space  weather  changes  such  as  occur  during  geomagnetic  storms. 


Fig.  14.  (T op)  GPS  ray  path  grids  for  plasma  tomography  using  ground  stations  (left)  and 

FedSat  (right). 

(Middle)  Structure  of  plasmasphere  during  magnetically  quiet  conditions. 
(Lower)  Structure  of  plasmasphere  during  magnetically  disturbed  conditions. 

The  white  lines  represent  the  Earth ’s  magnetic  field  lines. 

Fig.  14  shows  initial  results  from  an  ongoing  study  of  the  plasmasphere  using  FedSat  [31]. 
The  top  panel  compares  the  typical  ray  grid  obtained  during  a  FedSat  pass  with  that  obtained 
using  5  ground  stations  over  Australia.  It  is  evident  that  the  FedSat  observations  provide  a 
denser  grid  giving  more  uniform  coverage  of  the  region  being  investigated.  The  other  two 
panels  show  comparisons  between  FedSat  observations  of  the  plasmasphere  and  model 
predictions  for  the  same  solar  geophysical  conditions.  The  model  developed  by  Webb  [32], 
provides  a  good  description  of  our  current  knowledge  of  the  detailed  structure  and  behaviour 
of  the  plasmasphere.  The  middle  panel  shows  for  relatively  undisturbed  geophysical  and  solar 
conditions,  there  is  quite  good  agreement  between  model  and  observation  in  the  background 


structure  of  the  plasmasphere,  but  the  observations  show  a  trough  feature  not  predicted  by  the 
model. 

The  lower  panel  show  a  comparison  under  very  disturbed  conditions  and  the  observations 
show  quite  complicated  structures  not  predicted  by  the  model.  I  the  plasmasphere  the  plasma 
is  dominated  magnetic  field  effects  and  so  structures  can  be  expected  to  be  field-aligned. 
There  is  evidence  of  this  is  the  observations  between  1500  and  4000  km  altitude  and 
50°-70°  S.  However  the  plasmasphere  is  dominated  by  generally  reduced  plasma  density  and 
an  enhancement  that  stretches  upwards  across  the  field  lines.  This  suggests  that  electric  fields 
are  playing  a  significant  role  in  transporting  plasma  across  the  field  lines  by  E  x  B  drift. 

These  two  examples  illustrate  our  lack  of  detailed  knowledge  of  the  structure  and  dynamics  of 
the  plasmasphere  under  changing  space  weather  conditions  and  our  FedSat  studies  have  the 
potential  to  extend  significantly  our  understanding  of  this  region. 


Other  Research  Developments 

Significant  research  programs  are  also  being  undertaken  using  vertical  incident  radars 
(ionosondes)  in  Antarctica  and  at  La  Trobe’s  Bundoora  campus.  The  Antarctic  program  is  a 
collaborative  program  with  Dr  Ray  Morris,  Australian  Antarctic  Division  and  has  included 
detailed  mapping  of  ionospheric  convection  at  Casey,  Antarctica,  a  polar  cap  station  [34], 
[32]. 

The  experimental  radio  studies  are  supported  by  ray  tracing  studies  [35],  the  development  of 
new  inversion  techniques,  including  on  of  FedSat  occultation  observations  to  determine 
ionospheric  structure  [36]. 

A  new  initiative  being  led  by  Dr  Conde  is  the  planned  development  of  rocket  payloads  to 
measure  properties  of  the  atmosphere  up  to  100  km  using  sounding  rockets  being  developed 
by  the  University  of  Queensland  [37], 


Conclusions 

The  complexity  of  the  coupling  from  the  solar  wind  to  Earth’s  magnetosphere,  ionosphere  and 
thermosphere  results  in  a  great  variety  of  space  weather  phenomena  that  cannot  all  be  studied 
by  any  one  technique.  Progress  is  dependent  on  coordinated  studies  using  world-wide 
networks  of  satellite  and  ground  based  techniques  involving  different  instruments  and 
international  collaboration  between  various  scientific  groups.  In  this  paper  three  techniques, 
HF  ionospheric  radar,  large  Fabry-Perot  Interferometers,  and  satellite-to-satellite  tomography, 
used  by  the  Space  Physics  research  Group  at  La  Trobe  University,  have  been  described.  Some 
recent  results  have  been  described  and  the  leading  research  questions  which  each  technique 
can  address  have  been  presented. 
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This  paper  presents  a  backscatter  inversion  technique 
that  uses  the  elevation  angle  and  group  path  of  the 
received  signals. 
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More  Realistic  Profiles  can  be  obtained  by  fitting 
more  QP  segments: 


Chen,  J.,  J.  A.  Bennett  and  P.  L.  Dyson. 

Automatic  fitting  of  quasi-parabolic  segments  to  ionospheric  profiles  with  application 
to  ground  range  estimation  for  single-station  location, 

J.  Atmos.  Terr.  Phys.,  52,  277-288,  1990. 
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Relationship  between  Group  Path  and  Elevation  Angle  on  a  single  frequency 


Each  layer  is  described  by: 

Hence  inversion  requires  at  least  3 

•  Critical  frequency  points  per  layer 

•  Peak  Height 

•  Semi-thickness 


DctaJC 


*  friR  in^^ClCU^riaL£llL£Li^CLLLL^  (wriSfS  J1  ^ 

datfi  [BS§ 


rom  each  of  the  layer  echo  traces, 
roup  paths  P'1,P'2,.--,P'n  are  chosen 
ing  to  the  elevations 
espectively. 

ixist  explicit  equations  for  the  Group  path 
^xplessed  fri  terms  of  the  layer  parameters  and  the 
initial  elevation  angle. 

B  The  method  then  sets  out  to  find  a  set  of  values  foe 
the  layer  parameters,  which  yield  to  within  a 
specified  accuracy,  to  the  data  points  chosen. 


4 


5 


Table  2  Shows  the  convergence  of  the  E  layer  parameters  to  the  correct  result  after  each  iteration 


E  Layer  Iteration  No. 
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hmE  [km] 

91.00000 

101.1981 

100.9853 

101.0000 

101.0000 
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Synthesized  TIGER  Sea  echoes  during  Spring  of  2000 
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Synthesized  TIGER  Sea  echoes  during  Spring  of  2000 
Divergent  Power  Loss  along  Beam  4  for  12.0  MHz 
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ABSTRACT 

High  Frequency  (HF)  radar  land  and  sea  backscatter  inversion  to  obtain  ionospheric  electron  density  profiles  has  the 
potential  to  provide  remote  sensing  of  the  ionosphere  up  to  thousands  of  kilometers  from  the  transmitter/receiver 
location.  The  inversion  technique  requires  HF  radar  data  in  the  form  of  time  delay  or  group  path,  Pf ,  as  a  function  of 
elevation  angle,  p ,  from  each  of  the  ionospheric  layer  echo  traces.  This  HF  radar  data  is  inverted  to  obtain  a  multiple 
quasi-parabolic  segment  (QPS)  plasma  frequency  profile  of  the  true  ionosphere.  The  technique  can  be  used  with  both 
fixed  frequency  and  frequency  scanning  radars.  The  advantage  of  using  frequency  scanning  radar  is  that  the  down-range 
gradients  in  electron  density  can  be  determined,  and  spatial  variations  of  the  ionosphere  within  the  radar  coverage  area 
can  be  modeled  and  mapped.  In  this  paper,  a  prototype  HF  radar  backscatter  inversion  technique  is  exercised  using  real 
ionospheric  data  recorded  with  the  Tasman  International  Geospace  Environment  Radar  (TIGER). 


INTRODUCTION 

The  inversion  technique  developed  by  Norman  and  Dyson  (2005)  [1]  is  applied  to  real  backscatter  data  from  the 
Tasman  International  Geospace  Environment  Radar  (TIGER)  [2].  TIGER  is  an  Over-the-Horizon  Radar  (OTHR) 
situated  on  Bruny  Island  (43.4°S,  147. 2°E;  55°A)  and  transmits  southward  toward  the  austral  auroral  oval.  TIGER  was 
designed  to  study  the  physical  processes  occurring  in  geospace  including  the  aurora  and  magnetospheric  convection. 
TIGER  is  part  of  the  Super  Dual  Radar  Network  [3],  [4]  consisting  of  7  radars  in  the  Southern  Hemisphere  and  10  in 
the  Northern  Hemisphere. 

TIGER  transmits  HF  signals  where  some  of  the  signal  is  reflected  from  the  ionosphere  to  a  distant  location  on  the 
Earth’s  surface,  which  can  be  up  to  a  few  thousand  kilometers  from  the  transmitter/receiver  location.  Some  of  the  signal 
is  reflected  back  to  the  transmitter  location  along  the  same  path  via  reflection  from  the  ionosphere.  The  backscattered 
signals  contain  useful  information  regarding  the  state  of  the  intervening  ionosphere  at  the  time,  and  over  the  range  of 
the  returned  signal.  Especially  strong  backscatter  from  the  sea  is  recorded  during  the  day  time,  and  the  elevation  angles 
of  the  received  signals  are  determined  at  every  group  range  using  an  interferometer  sub-array  [5]. 

Norman  and  Dyson  [1]  describe  the  backscatter  inversion  technique  in  detail  where  the  theory  and  accuracy  of  the 
technique  are  shown.  Their  inversion  technique  requires  HF  radar  data  in  the  form  of  time  delay  or  group  path,  P' ,  as  a 
function  of  elevation  angle,  p ,  from  each  of  the  ionospheric  layer  echo  traces.  This  inversion  technique  determines 
quasi-parabolic  ionospheric  layer  parameters  [6]  for  each  of  the  layer  echo  traces  presented  in  the  form  of  profiles  of 
group  path  versus  elevation  angle.  At  least  3  group  path  -  elevation  angle  data  pairs  are  required  from  each  layer  echo 
trace,  and  in  general  the  more  data  points  the  better  the  inversion  results  will  be.  The  maximum  elevation  of  each  of  the 
ionospheric  layer  echo  traces  is  also  a  required  input,  and  is  used  to  seed  the  first  guess  of  the  layer  parameters.  This 
data  is  then  inverted  to  form  a  multiple  quasi-parabolic  segment  (QPS)  model  ionospheric  profile.  The  QPS  model 
requires  5  QP  segments  to  represent  the  E,  FI,  and  F2  ionospheric  layers. 
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TIGER  has  a  main  array  of  16  log-periodic  antennas  and  the  use  of  a  phasing  matrix  permits  the  main  beam  to  be 
steered  in  16  different  directions  separated  by  3.24°  to  achieve  an  azimuthal  coverage  of  -52°.  TIGER  is  capable  of 
illuminating  the  ionosphere  over  a  wide  range  of  elevation  angles  and  azimuth  angles.  TIGER  is  frequency  agile  and  in 
the  most  common  modes  of  operation  automatically  selects  a  frequency  appropriate  for  observations.  Once  a  frequency 
is  chosen,  the  radar  operates  at  that  frequency  for  extended  periods  of  time.  Thus  the  inversion  technique  has  been 
developed  for  a  fixed  operating  frequency.  However,  by  varying  the  operating  frequency  of  the  HF  radar  it  is  possible 
to  increase  the  ionospheric  coverage  area  for  the  inversion  technique  and  determine  the  down-range  gradients  of  plasma 
frequency.  Thus  the  inversion  technique  can  be  performed  on  multi-frequency  data  collected  on  each  of  the  16  beams  to 
construct  realistic  3-D  mappings  of  ionospheric  electron  density. 

Norman  and  Dyson  [1]  showed  an  example  of  the  technique  using  real  backscatter  data  consisting  of  three  ionospheric 
layer  echo  traces.  Their  example  used  group  path  versus  elevation  angle  profiles  recorded  at  a  single  operating 
frequency.  In  this  paper,  the  example  chosen  consists  of  group  path  versus  elevation  angle  profiles  recorded  over  a 
range  of  operating  frequencies,  and  thus  enabling  a  description  of  horizontal  gradients  in  the  ionosphere  to  be 
determined,  at  least  immediately  downrange  from  the  radar  site.  In  this  initial  example,  only  data  recorded  on  a  single 
beam  under  relatively  simple  conditions  will  be  presented  to  illustrate  the  capabilities  of  the  technique. 


THEORY 

The  backscatter  inversion  technique  is  described  in  detail  in  [1]  where  a  QPS  model  is  used  to  describe  the  plasma 
frequency  (proportional  to  the  square  root  of  electron  density)  as  a  function  of  the  height  of  the  ionosphere  near  the 
midpoint  of  the  propagation  path.  The  model  is  constructed  from  the  backscattered  signal  from  each  of  the  layer  echo 
traces,  which  are  in  the  form  of  group  path  versus  elevation  angle  profiles. 

Up  to  5  QP  segments  are  required  to  represent  a  typical  daytime  ionosphere.  Three  of  the  QP  segments  represent  the  3 
ionospheric  layers,  namely  the  E,  FI  and  F2  layers,  and  2  QP  segments  are  joining  segments  between  the  layers  and 
smoothly  attach  to  the  layer  QP  segments.  A  QP  layer  consists  of  only  3  layer  parameters,  namely  the  peak  plasma 
frequency  of  the  layer,  the  height  of  maximum  electron  density  of  the  layer,  and  the  semi-thickness  of  the  layer. 

In  this  paper,  the  inversion  technique  is  used  on  a  backscatter  ionogram  example  containing  a  section  of  echo  trace 
identified  as  single  F-region  layer  echo  trace.  Thus  the  inversion  technique  will  produce  results  showing  a  single  QP 
segment.  The  backscatter  ionogram  provides  the  group  path  versus  elevation  angle  for  a  range  of  operating  frequencies. 
By  applying  the  inversion  technique  on  the  group  path  versus  elevation  angles  at  each  operating  frequency,  the  down- 
range  ionospheric  profile  can  be  modeled  and  mapped. 

The  inversion  technique  requires  at  least  three  data  points  in  the  form  of  group  path  versus  elevation  angle  from  each  of 
the  layer  echo  traces.  The  maximum  elevation  angle  of  each  of  the  layer  echo  traces,  /?max ,  is  a  required  input  and  is 

used  to  seed  the  initial  guess  ionospheric  ray  parameters  (fc,rb,rm).  The  seed  values  of  rm  and  rb  can  be 
incremented  so  that  the  process  homes  in  to  the  best  possible  solution. 

The  seed  layer  parameters  are  then  substituted  into  analytic  expressions  for  the  group  path  where  the  computed  group 
paths,  P'cl ,  P'2 ,. . Pcn  will  most  likely  differ  from  the  real  values  P{ ,  P2  ,. . .,  P^ .  Let  the  corresponding  differences 
between  the  real  and  computed  values  be  A P( ,  A P2,...,  A P„  ,  respectively.  Then,  to  a  first  approximation  the  amount 
Sxt  =  ( A fc  ,  A  rb  ,  A  rm  )  by  which  the  assumed  layer  parameters  should  be  incremented,  so  that  A P{ ,  AP2' AP'  are 
minimised,  is  given  by 


A P\f)  =  I.K^-dXi 

i  OX: 


Using  the  same  notation  as  [5]  the  equation  can  be  re-written  in  matrix  form: 
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or  more  simply  as  Z  =  ST 
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is  a  (1,  n)  matrix, 
is  a  (3,  n)  matrix, 
is  a  (3,  1)  matrix. 


The  inversion  of  this  equation  is  T  =  [STS]  1 STZ 


It  should  be  noted  that  these  equations  may  be  found  in  [1]  where  ST  is  the  transpose  matrix  of  S  and  the  matrix  ST  S 
is  a  square  matrix  capable  of  inversion. 

The  assumed  layer  parameter  values  are  then  incremented  by  A fc ,  A rb  and  A rm  .  Thus  a  new  set  of  layer  parameters  is 
determined.  The  process  is  repeated  until  a  set  of  values  for  the  layer  parameters  produces  group  path  values,  which 
yield  to  within  a  specified  accuracy,  to  the  chosen  real  group  path  versus  elevation  angle  data  points. 

Using  multiple  operating  frequencies  the  inversion  technique  can  be  applied  to  data  points  received  from  each 
frequency  highlighting  different  regions  of  the  ionosphere.  An  overall  mapping  of  the  ionosphere  in  the  coverage  area 
of  the  radar  can  then  be  performed. 


TIGER  RESULTS 

Fig.  1  is  a  TIGER  backscatter  ionogram  recorded  on  the  magnetic  meridian  pointing  beam  4  commencing  at  21:00  UT 
on  20th  February  2000.  The  backscatter  ionogram  displays  the  signal-to-noise  ratio  (SNR)  of  the  received  signal  in  dB 
versus  transmitter  frequency  and  group  range.  The  group  range,  or  group  path,  is  obtained  from  the  time  delay  of  the 
echo,  converted  to  a  distance  assuming  propagation  at  the  speed  of  light  in  a  vacuum.  TIGER  was  scanning  in 
frequency  from  8  to  20  MHz  in  steps  of  approximately  0.1  MHz.  The  integration  time  was  3  s  per  frequency;  hence  the 
backscatter  ionogram  took  6  minutes  to  acquire.  The  SNR  values  were  recorded  at  steps  of  45  km  using  Gaussian-like 
pulses  of  width  300  ps,  equivalent  to  45  km  in  range. 

SuperDARN  radars  also  measure  the  line-of-sight  Doppler  velocity,  spectral  widths,  and  elevation  angles  of  echoes  [4]. 
Figs.  2  to  4  are  the  corresponding  backscatter  ionograms  showing  these  parameters.  A  real  time  radar  algorithm 
automatically  tags  echoes  as  emanating  from  the  ground  or  sea  if  they  have  small  Doppler  speed  and  spectral  width, 
within  certain  error  limits.  The  main  echo  trace  which  dominates  the  ionograms  has  Doppler  speeds  and  spectral  widths 
less  than  ~15  m  s'1.,  and  thus  represents  1-hop  F-region  backscatter  from  the  Southern  Ocean.  On  the  other  hand,  the 
echo  trace  with  high  Doppler  speed  from  12-13  MHz  and  2300-2800  km  represents  1.5-hop  backscatter  from  a  patch 
of  ionospheric  irregularities. 
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Fig.  1.  OTHR  ionogram  showing  the  backscatter  power  (SNR)  recorded  on  TIGER  beam  4  at  21:00  UT  (07:00  LT)  on 
the  20th  February  2000. 
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Fig.  2.  The  same  as  Fig.  1  except  the  backscatter  ionogram  shows  the  FOS  Doppler  velocity. 
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.  The  same  as  Fig.  1  except  the  backscatter  ionogram  shows  the  Doppler  spectral  widths. 
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Fig.  4.  The  same  as  Fig.  1  except  the  backscatter  ionogram  shows  the  elevation  angles. 


Both  Figures  1  and  4  show  that  there  is  complicated  structure  in  the  echo  traces  beyond  2000  km  in  range.  There  is  an 
echo  trace  that  starts  at  about  12  MHz  with  a  leading  edge  at  about  2200  km  range.  At  higher  frequencies  the  signals  are 
stronger  and  this  trace  merges  with  the  main  trace  whose  leading  edge  increases  linearly  with  frequency.  This  echo 
structure  suggests  that  there  was  a  dramatic  change  in  the  structure  of  the  ionosphere  at  around  2000  km  distant  from 
the  radar.  The  form  of  the  inversion  technique  presented  here  assumes  that  the  ray  paths  associated  with  sea  scatter  are 
symmetric  to  first  order,  so  it  will  not  give  a  good  representation  of  sharp,  almost  discontinuous  gradients,  which  may 
be  associated  with  the  occurrence  of  the  additional  echo  trace  beyond  2000  km  range.  A  more  sophisticated  inversion 
method  that  could  deal  with  such  structures  was  described  by  Norman  and  Dyson  (2005)  which  requires  the  Segment 
Method  Analytic  Ray  Tracing  SMART  technique  [7  &  8].  Consequently  for  our  purpose  here  of  demonstrating  the 
basic  inversion  technique,  the  analysis  has  been  limited  to  the  main  F  layer  trace  up  to  a  frequency  of  14  MHz. 


INVERSION  RESULTS 

Fig.  5  represents  the  resultant  structure  of  the  ionosphere  along  the  magnetic  meridian  pointing  beam  4  when  the  group 
path  versus  elevation  angle  profiles  contained  within  the  main  F  layer  trace  in  the  backscatter  ionogram,  Fig.  4,  were 
fed  into  the  inversion  program.  The  plot  shows  the  plasma  frequency  in  kHz  versus  the  ground  range  and  true  height. 
The  exact  value  of  the  height  of  the  peak  electron  density,  hm F2,  has  been  highlighted  using  bold  horizontal  bars,  and 
the  numbers  located  just  beneath  the  base  of  the  F  layer  represent  every  fourth  value  of  the  operating  frequency  in  kHz. 
A  simple  Chapman  layer  was  included  to  describe  the  ionosphere  above  hm F2. 

The  Fig.  5  results  have  been  pixilated  using  a  nominal  resolution  of  1  km  in  true  height,  and  the  width  in  ground  range 
is  variable  but  controlled  by  the  step  size  used  in  the  frequency  scan.  The  inversion  applied  to  a  single  frequency  yields 
a  vertical  profile  of  the  ionosphere  nominally  located  near  to  the  mid-point  of  refraction,  but  it  is  actually  related  to  a 
significant  spread  of  range  determined  by  oblique  propagation  through  the  various  layers  of  the  ionosphere.  We  are  in 
the  process  of  quantifying  this  effect  which  will  depend  on  the  details  of  the  inversion  procedure. 

In  the  absence  of  low  frequency  measurements  providing  group  range  versus  elevation  angle  data  points  for  the  E  layer, 
a  simple  model  E  layer  was  included  in  the  inversion.  The  chosen  parameters  were/0E=2.5  MHz,  hm E=115  km,  and  a 
semi- thickness,  ymE=15  km.  The  E-layer  parameters  were  not  varied  with  ground  range,  but  in  principle  they  could 
have  been  by  obtaining  values  from  a  standard  ionospheric  model.  The  main  effect  of  including  an  E  layer  in  the 
inversion  was  to  reduce  hm F2  for  the  F  layer,  especially  at  the  furthest  ranges.  Increasing  the  peak  density  and  semi¬ 
thickness  of  the  E-layer  would  reduce  hm F2  further  than  shown  here. 

Fig.  5  shows  how  both  hm F2  and  the  peak  electron  density  of  the  F  layer,  f0 F2,  fluctuated  with  ground  range,  but  overall 
they  increased  further  away  from  the  radar.  At  ground  ranges  <400  km,  hm F2  was  -270  km  and/0F2  was  -6.3  MHz  and 
increasing.  At  the  furthest  ground  ranges  >1200  km,  hm F2  reached  -440  km  and/0F2  reached  >9  MHz.  The  ionosphere 
at  furthest  ranges  was  unusually  high  and  dense,  and  the  reality  of  down-range  gradient  needs  direct,  independent 
corroboration. 

The  Australian  Space  Weather  Agency  (http://www.ips.gov.au)  operates  ionosondes  at  Hobart  (42.9°S,  147. 3°E), 
located  close  to  the  TIGER  radar  site  (43.4°S,  147. 2°E),  and  at  Macquarie  Island  (54.5°S,  158. 9°E;  -65 °A),  located 
downrange  and  just  to  the  east  of  the  TIGER  radar  field  of  view.  The  geomagnetic  (geographic)  latitude  of  Macquarie 
Island  corresponds  to  a  ground  range  of  -1035  km  on  TIGER  beam  4,  but  is  separated  by  a  great  circle  distance  of 
-1010  km.  The  distribution  of  plasma  density  will  be  partly  controlled  by  electrodynamics  and  solar  illumination  at 
auroral  latitudes. 


The  Hobart  and  Macquarie  Island  ionosondes  recorded  hourly  values  of/0F2=6.3  MHz  and  7.3  MHz,  respectively,  at 
21:00  UT  on  20  February  2000  (see  http://www.ips.gov.au/World  Data  Centre/1/2).  In  Fig.  5,/0F2  was  -6.3  MHz  at 
the  closest  ground  ranges,  in  agreement  with  the  Hobart  ionosonde.  However,  fQ F2  varied  between  -8  to  9  MHz  at 
ground  ranges  1000  to  1300  km.  Hence  the  down-range  increase  of  plasma  density  was  greater  than  implied  by  the 
ionosonde  data.  However,  Macquarie  Island  was  located  far  away  from  beam  4,  and  large,  dynamic  variations  in  plasma 
density  are  expected  at  auroral  latitudes.  Although  the  IRI  model  is  not  accurate  at  these  latitudes,  it  indicates  that  hm E2 
should  have  steadily  increased  with  latitude,  but  only  from  -273  to  280  km  over  the  ground  range  shown  in  Fig.  5. 
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Fig.  5.  The  down-range  variation  of  ionospheric  plasma  frequency  determined  by  inverting  the  backscatter  ionogram 
data  recorded  on  TIGER  beam  4  at  21:00  UT  (07:00  LT)  on  the  20th  February  2000. 


No  E-layer  parameters  were  available  for  the  Hobart  and  Macquarie  Island  ionosondes  at  21:00  UT  on  20  February,  but 
on  the  preceding  day,  fQ E  was  2.65  MHz  and  2.85  MHz  at  the  two  stations,  respectively.  Including  a  uniform  model  E 
region  with/0E=2.5  MHz  in  the  inversion  was  conservative,  and  this  may  partly  account  for  the  large  values  of  hm F2 
and/0F2  obtained  at  the  furthest  ranges. 


TIGER  observes  a  remote  and  poorly  known  region  of  the  Earth’s  ionosphere,  so  development  of  this  technique  would 
provide  valuable  new,  routine  information  on  the  electron  density  structure.  Fig  5  shows  an  increase  in  foF2  with 
latitude,  i.e.  distance  from  the  radar.  We  note  that  Horvath  and  Essex  [9]  also  found  the  electron  density  to  increase  with 
latitude  in  this  region.  They  analyzed  total  electron  content  (TEC)  measurements  and  discovered  a  dramatic  “build  up  of 
ionization”  located  just  equator  ward  of  the  daytime  main  ionospheric  trough. 

The  inversion  algorithm  can  quickly  home  in  on  layer  parameters  with  a  precision  far  exceeding  that  obtainable  from 
observations  [1].  The  errors  in  the  inversion  are  probably  dominated  by  errors  in  the  measured  group  range  versus 
elevation  angle  profiles.  Only  the  4  group  path  versus  elevation  data  points  with  maximum  SNR  were  automatically 
chosen  at  each  operating  frequency  for  the  present  analysis.  Random  errors  in  the  elevation  angle  data,  combined  with 
the  effect  of  the  criteria  used  to  automatically  select  the  data  points,  probably  dominate  the  fluctuations  shown  in  Fig.  5. 


SUMMARY  AND  FUTURE  WORK 

This  new  backscatter  inversion  technique  offers  a  means  of  sensing  remote  areas  of  the  ionosphere  where  it  is  difficult 
to  obtain  vertical  ionograms.  The  down-range  gradients  in  electron  density  in  the  coverage  area  of  a  radar  network, 
where  the  group  range  versus  elevation  angle  profiles  are  recorded,  can  be  readily  determined  using  this  backscatter 
inversion  technique.  The  procedure  can  be  fully  automated,  as  it  was  in  the  present  example,  and  the  resultant 
ionospheric  layer  parameters  and  ray  tracing  results  were  obtained  in  a  matter  of  seconds.  Applying  the  technique  to 
multi- frequency,  multi-beam  radar  measurements,  3-D  maps  of  the  ionosphere  within  the  radar  coverage  area  can  be 
obtained. 
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[i]  Methods  of  inverting  high-frequency  radar  land  and  sea  backscatter  to  obtain 
ionospheric  vertical  profiles  offer  an  important  means  of  remote  sensing  the  ionosphere  up 
to  thousands  of  kilometers  from  the  transmitter/receiver  location.  This  new  inversion 
technique  requires  the  received  elevation  angle  as  well  as  the  time  delay  or  group  path 
of  the  backscattered  signal.  The  technique  can  be  used  with  both  fixed- frequency  4nd 
swept-frequency  radars.  The  advantage  of  using  swept-frequency  radar  is  that  the 
downrange  gradients  in  electron  density  can  more  readily  be  determined. 
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1.  Introduction 

[2]  The  case  we  consider  is  ground  backscatter  signals 
received  from  distant  locations  on  the  Earth’s  surface  via 
reflection  by  the  ionosphere.  Both  the  transmitted  and 
received  signals  traverse  the  ionosphere  and  these  signals 
contain  useful  information  regarding  the  state  of  the 
intervening  ionosphere  at  the  time,  and  over  the  range 
of  the  returned  signal  which  can  be  up  to  a  few  thousand 
kilometers  from  the  transmitter/receiver  location,  ifiphis 
paper  a  new  HF  radar  backscatter  inversion  method  is 
presented.  It  is  similar  in  concept  to  that  of  Rao  [1910, 
Bertel  et  al.  [1987],  and  Norman  [2003],  where  inversion 
techniques  were  derived  for  backscatter  ionogram  inver¬ 
sion.  This  new  inversion  technique  requires  HF  radar 
data  in  the  form  of  time  delay  or  group  path,  P',  as  a 
function  of  elevation  angle,  (3,  hop  each  of  the  iono¬ 
spheric  layer  echo  traces.  This  data  is  then  inverted  to 
form  a  multiple  quasi-parabolic  segment  ionospheric 
profile,  using  the  quasi-parabolic  segment  (QPS)  model 
derived  by  Dyson  and  Bennett  [19l8].  This  QPS  model 
uses  five  QPSs  to  represent  the  E ,  Fx  and  F2  ionospheric 
layers. 

[3]  The  characteristics  of  the  backscattered  echoes 
depend  on  the  radar  system  characteristics  and  iono¬ 
spheric  characteristics  as  well  as  the  backscattering 
properties  of  the  ground.  Important  radar  characteristics 
are  the  wavelength,  pulse  length,  antenna  beam  width, 
antenna  gain,  and  transmitted  power.  HF  radars  that 
utilize  ionospheric  propagation  at  oblique  incidence  are 
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usually  capable  of  illuminating  the  ionosphere  over  a  45 
relatively  wide  range  of  elevation  angles,  giving  rise  to  46 
the  possibnw  of  ground  backscatter  being  detected  over  47 
iyyide  span  of  ranges.  The  radar  characteristics  largely  48 
affect  the  signal  strength  of  echoes  and  the  ionospheric  49 
structure  determines  the  Pr  elevation  angle  characteristics  50 
at  each  frequency  so  throughout  this  investigation  we  51 
^rasider  only  these  ionospheric  effects.  Displays  of  52 
echoes  detected  by  ionospheric  HF  radars  are  commonly  53 
called  ionograms  and  here  we  are  primarily  concerned  54 
with  fixed-frequency  backscatter  ionograms  displaying  55 
echo  traces  formed  by  the  variation  of  P'  with  elevation  56 
angle.  57 

[4]  The  inversion  technique  was  tested  using  real  58 

backscatter  data  from  the  Tasman  International  Geospace  59 
Environment  Radar  (TIGER)  [Dyson  and  Devlin,  2000].  60 
TIGER  is  an  oblique  incidence,  ionospheric  backscatter  61 
radar,  designed  to  study  the  physical  processes  occurring  62 
in  space  weather  systems  such  as  radio  aurora.  It  is  63 
situated  on  Bruny  Island  (43.4°S,  147.2°E;  55°A)  and  64 
transmits  southward  and  is  part  of  the  Super  Dual  Radar  65 
Network  [Greenwald  et  al,  1985,  1995]  consisting  of  six  66 
radars  in  the  Southern  Hemisphere  and  nine  in  the  67 
Northern  Hemisphere.  Throughout  this  investigation  we  68 
concentrate  on  the  ionospheric  effects  on  the  propaga-  69 
tion;  that  is,  it  is  assumed  that  the  radar  characteristics  are  70 
known  and  can  be  accounted  for  separately.  71 

[5]  This  new  inversion  technique  derives  quasi-para-  72 
bolic  ionospheric  layer  parameters  from  the  layer  echo  73 
traces  in  the  form  of  group  path  and  elevation  angle.  At  74 
least  three  group  path  elevation  angle  data  points  are  75 
required  from  each  layer  echo  trace.  The  maximum  76 
elevation  of  each  of  the  ionospheric  layer  echo  traces  77 
is  also  a  required  input.  The  inversion  technique  is  based  78 
on  using  a  fixed  operating  frequency.  However,  by  79 
varying  the  operating  frequency  of  the  HF  radar  it  is  80 
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possible  to  determine  the  downrange  ionospheric  gra¬ 
dients  and  thus  produce  more  realistic  ionospheric 
mappings. 

2.  Quasi-Parabolic  Segment  QPS  Model 

[6]  The  QP  layer  defined  by  Croft  and  Hoogasian 
[1968]  is  given  by 


Ne 


fr  —  rm\  2  /rb\2 

Nm 

i-  — -  - 

< 

V  ym  J  \rJ 

.  0, 

rb  <  r  <rm 


otherwise. 


NJ  =  *e~ 


and  the  F  layer  is 

Np  =  ap  — 

where  a  =  Nm  and  b  =  Nm  ( rb/y m)2  and 

bj  =  —rFbF  ^1  — — 

where  the  joining  point,  rc,  at  which  the  joining  layer 
connect  to  the  F  layer  may  be  written  as 


Wii1) 

aF  ~  aE  +  bp  —  1^ 


[9]  If  an  Fi  layer  is  present  then  QPS  joining  layers  are  119 

used  between  the  Fx  and  F2  QPS  layers  and  between  the  120 
Fi  and  E  QPS  layers.  121 

[10]  The  work  of  Dyson  and  Bennett  [1988]  provides  122 

exact  expressions  for  calculating  ray  parameters  such  as  123 
the  group  path,  P',  for  propagation  in  a  spherically  124 
stratified  ionosphere  consisting  of  a  single  QP  layer.  125 
The  equation  for  calculating  the  group  path  may  be  126 
written  as  L  >  127 


i  r*  X. 

<N 

II 

1  n  -  ym) 

r 


rdr 


r\fr 


V 'Ar2  +  Br  +  C 


(1) 


where 

Ne  electron  density  at  a  radial  distance  r  from  the 
Earth’s  center; 

Nm  maximum  electron  density  at  the  radial  height  rm; 

rb  radial  base  height  of  the  ionospheric  layer; 

ym  layer  semilayer  thickness. 

[7]  The  QPS  model,  presented  in  the  work  of  Dyson 
and  Bennett  [1988],  consists  of  smoothly  joine<%QP 
segments,  where  QPSs  are  used  to  describe  each  of  ti 
ionospheric  layers  at  and  immediately  below  their  peaks 
and  additional  QPSs  are  used  as  joining  segments  that 
smoothly  join  together  the  QP  descriptions  of  the  layer 
peaks.  The  equations  describing  the  QPjg  ionospheric 
model  may  be  written  as  follows. 

[8]  The  E  layer  is 

Ne  =  ap  —  bp(^  1  - 
the  joining  layer  is 


where  rt  is  the  radial  height,  from  the  Earth’s  center,  at  129 
which  the  raypsmflected  and  130 

^p2  —  r2  cos2  (3  =  Ar 2  +  Br  +  C, 

where  p  is  the  refractive  index  and  (3  is  the  elevation  132 
angle.  133 

[11]  Dyson  and  Bennett  [1988]  presented  explicit  equa-  134 
tions  for  the  ray  parameters  in  their  QPS  model  for  a  135 
spherically  stratified  ionosphere,  where  if  the  ray  prop-  136 
agates  through  n  segments  the  group  path  may  be  written  137 

as  138 


P'  =  2 


-  r20  cos2  (3  —  rQ  sin  (3  +  ^  (Un  -  Ln 


where  Un  and  Ln  represent  the  values  of  the  integral  in  140 
equation  (1)  at  the  upper  and  lower  bounds.  141 

3.  New  Radar  Backscatter  Inversion  142 

Technique  143 

[12]  The  method  of  inversion  of  the  backscattered  echo  144 

traces  presented  here  requires  the  determination  of  the  145 
three  layer  parameters  for  each  of  the  quasi-parabolic  146 
layers,  namely  the  critical  frequency  fc,  base  height  rb,  147 
and  height  of  maximum  electron  density  rm.  148 

[13]  The  inversion  technique  begins  with  inversion  of  149 

the  backscattered  signal  refracted  from  the  E  layer.  Once  150 
the  parameters  for  the  E  layer  have  been  calculated  the  15 1 
inversion  technique  inverts  the  Fi  layer  echo  trace  using  152 
the  E  layer  parameter  results  as  well  as  the  QPS  equa-  153 
tions  for  the  joining  layer,  which  smoothly  joins  the  E  154 
layer  and  the  Fx  layer.  The  process  is  then  repeated  for  155 
the  F2  layer  where  the  layer  parameters  already  deter-  156 
mined  for  the  E  and  Fx  layers,  as  well  as  the  equations  157 
for  the  joining  layer  that  joins  the  Fx  and  F2  layers  are  158 
required  for  the  inversion  process.  159 

[14]  The  inversion  technique  requires  n  (where  n  is  at  160 
least  3)  data  points  from  each  of  the  layer  echo  traces,  161 
where  the  n  group  paths  chosen  P[,  Pf. .  .,Pn  correspond  162 
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Elevation  angle  [deg] 

Synthesized  group  path  versus  elevation  plot  using  the  three-layer  i 

© 


Figure  1. 

with  the  layer  parameters  given  in  the  text. 


ionospheric  model 


to  the  n  elevations  (3i,  (32  -  -  .(3„,  respectively.  The  method 
then  sets  out  to  find  a  set  of  values  for  the  layer 
parameters,  which  yield  to  within  a  specified  accuracy, 
to  the  data  points  chosen.  To  accomplish  this  an  initial 
guess  of  the  ray  parameters  (rc,  rb,  rm)  is  made. 

[15]  The  maximum  elevation  angle  of  each  of  the  layer 
echo  traces,  (3max  is  a  required  input  and  is  used  to 
determine  the  initial  guess  ionospheric  parameters  and 
aids  in  stabilizing  the  technique.  A  raypath  having  an 
elevation  angle  of  (3max  reaches  an  apogee  height  equal  to 
the  peak  height,  or  maximum  height,  rm  of  the  layer 
Raypaths  with  elevation  angles  greater  than  (3max  peng- 
trate  the  ionospheric  layer. 

[16]  Then,  using  the  maximum  elevation  angle  the 
layer  echo  trace  and  an  educated  guess  ofihe  peak  height 


Table  1.  Three  Data  Points  Chose: 
Echo  Traces 


Data  Points 

Elevation  deg 

Group  Path  PE,  km 

1 

EUyer 

7.0 

9.0 

1361.49009944569 

2 

1168.63887920563 

3 

1033.35342695743 

PmaxZ? 

12.0 

1 

F]  layer 

19.0 

1404.10692246982 

2 

19.2 

1348.36189670986 

3 

19.4 

1319.28979088134 

PmaxFl 

19.7 

1 

F2  layer 

22.0 

1238.53947948455 

2 

23.0 

1187.14934478616 

3 

24.0 

1155.70487464130 

PmaxF2 

29.3 

A 

■  base  hei 


rm  and  the  base  height  rb  of  the  ionospheric  layer,  the 
maximum  plasma  frequency  of  the  layer  can  be  deter¬ 
mined  using  the  equation 


fc 


-7T  COSZ 


178 

179 

180 


initial  guess  values  of  rm  and  rb  can  be  incremented  182 
so  that  the  process  homes  in  to  the  best  possible  solution.  183 

[17]  These  new  layer  parameters  are  then  substituted  184 
into  the  analytic  expressions  for  the  group  path  for  the  185 
quasi-parabolic  layer.  Let  the  computed  minimum  group  186 
paths,  which  most  likely  will  differ  from  the  real  values  187 
P\,  Pr2 ,  . . .,  P'n ,  be  P'ci ,  P'c 2,  . . P'cm  respectively,  and  the  188 
corresponding  differences  between  the  real  and  comput-  189 
ed  values  be  AP\,  A  P'2,  . . .,  A P'm  respectively.  Then,  to  a  190 
first  approximation  the  amount  =  (Afc,  Arb,  A rm)  by  191 
which  the  assumed  layer  parameters  should  be  incre-  192 
mented,  so  that  AP[,  A P2,  . . .,  AP'n  are  a  minimum  is  193 


Table  2.  Convergence  of  the  E  Layer  Parameters  to  the  t2.1 
Correct  Result  After  Each  Iteration 


E  Layer  Iteration 

t2.2 

t2.3 

1 

2 

3 

4 

5 

foE,  MHz 

3.253448 

3.344400 

3.319285 

3.319999 

3.320000 

t2.4 

ymE,  km 

13.00000 

14.53291 

13.98480 

13.99999 

14.00000 

t2.5 

hmE,  km 

91.00000 

101.1981 

100.9853 

101.0000 

101.0000 

t2.6 

P'ci,  km 

1255.018 

1359.362 

1361.448 

1361.490 

1361.490 

t2.7 

P'c 2,  km 

1070.978 

1166.943 

1168.593 

1168.639 

1168.639 

t2.8 

Pci,  km 

943.9176 

1031.923 

1033.308 

1033.353 

1033.353 

t2.9 

APj,  km 

106.4718 

2.12810 

4.188E-2a 

2.971E-5 

2.87E-9 

t2.10 

AF2,  km 

97.66112 

1.69550 

4.579E-2 

3.412E-5 

1.56E-10 

t2.11 

AP'3,  km 

89.43584 

1.43042 

4.553E-2 

3.631E-5 

-.39E-9 

t2.12 

aRead4.188E-2  as  4.188  x  10“ 


t2.13 
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t3.1  Table  3.  Convergence  of  the  F\  Layer  Parameters  to  the  Correct  Result  After  Each  Iteration 
t3.2  F\  Layer  Iteration 


t3.3 

1 

2 

3 

4 

5 

6 

t3.4 

foFl,  MHz 

4.716269 

4.736789 

4.744927 

4.749492 

4.749995 

4.750000 

t3.5 

ymFl,  km 

30.40000 

41.18461 

41.70168 

42.85851 

42.99864 

43.00000 

t3.6 

hmFl,  km 

152.0000 

163.7422 

163.7117 

163.9681 

163.9997 

164.0000 

t3.7 

P'd,  km 

1176.024 

1329.081 

1319.607 

1319.270 

•  1319.289 

1319.290 

t3.8 

P'd,  km 

1188.481 

1363.036 

1349.220 

1348.374 

1348.362 

1348.362 

t3.9 

Pfc 3,  km 

1214.872 

1437.361 

1408.569 

1404.416 

1404.109 

1404.107 

t3.10 

APfl,  km 

143.2662 

-9.791348 

-0.3176219 

1.9871E-2a 

3.1598E-4 

2.6038E-8 

t3.11 

AP2,  km 

159.8804 

-14.67364 

-0.8576807 

-1.1789E-2  ^ 

8.8950E-5 

-1.318E-9 

t3.12 

AP3,  km 

189.2347 

-33.25370 

-4.462045 

-0.3086568 

^k-2^437E-3 

-2.223E-7 

t3.13 


194 

195 


197 


198 


aRead  1.9871E-2  as  1.9871  x  10“2. 

given  by  (refer  also  to  Bertel  et  al.  [1987]  and  Norman 
[2003]) 

l  1 

and  may  be  rewritten  in  matrix  form  as 

fdP'l 


where 


201 


|AP'|  = 

or  in  expanded  form  as 


dx, 


\dxi\ 


A  P[ 


A P', 


A  P' 


dT\ 

drJp, 


dd\ 

dr  Jr 


dr_\ 

drJts„ 


or  more  simpr 
Z=  [AE*']  (l,n)  matrix; 

S  =  M;  (3,n)matrix; 

K  =  [<ix  1 )  matrix. 

The  inversion  of  this  equation  is  K  =  [Nr  S]~l  ST  210 
ST  is  the  transpose  matrix  of  S  and  the  matrix  211 
square  matrix  capable  of  inversion.  212 

The  assumed  layer  parameter  values  are  then  213 
remented  by  A fc,  A rb  and  A rm.  The  entire  procedure  214 
begins  again  with  the  new  layer  parameter  values  until  215 
the  differences  in  group  path  APi,  AP/2,  . .  .Ai^  converge  216 
to  a  small  specified  minimum,  thereby  obtaining  the  final  217 
solution  of  the  layer  parameters  for  that  ionospheric  218 
layer.  219 

[20]  Once  the  layer  parameters  have  been  calculated  220 
the  next  ionospheric  layer  parameters  are  then  solved  for  221 
using  the  technique  above  and  the  layer  parameters  222 
already  determined.  For  example  let  us  assume,  that  223 
the  E  layer  parameters  rbE,  rE,  foE  have  already  been  224 
evaluated  using  the  method  shown  above.  The  peak  225 
elevation  angle  of  the  Fx  layer  echo  trace,  (3maXjpi,  may  226 


t4.1 

t4.2 

t4.3 

Table  4. 

Convergence  of  the  F2  Lay 

er  Parameters  to  the  Correct  Result  After  Each  Iteration 

V 

F2  Layer  Iteration 

1 

2 

3 

4 

5 

t4.4 

foF2,  MHz 

6.450973 

6.450973 

6.445580 

6.449979 

6.450000 

t4.5 

ymF2,  km 

54.50000 

72.66666 

67.77551 

67.99893 

68.00000 

t4.6 

hmF2,  km 

218.0000 

218.0000 

213.8218 

213.9990 

214.0000 

t4.7 

P'd,  km 

1414.500 

1259.923 

1238.608 

1238.538 

1238.539 

t4.8 

Ki,  km 

1337.601 

1208.438 

1187.177 

1187.148 

1187.149 

t4.9 

P'c3,  km 

1277.102 

1178.111 

1155.730 

1155.704 

1155.705 

t4.10 

AP[,  km 

-175.9608 

-21.38345 

-6.88078E-2a 

1.278789E-3 

2.8484237E-8 

t4.11 

AP2,  km 

-150.4515 

-21.28879 

-2.71656E-2 

1.119846E-3 

2.0539801E-8 

t4.12 

AP3,  km 

-121.3973 

-22.40620 

-2.52626E-2 

9.308 160E-4 

1.7825414E-8 

t4.13 

“Read  6.88078E-2  as  6.8807  x  10“2. 
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be  determined  directly  from  the  backscattered  trace  and 
at  least  three  data  points  are  required  from  the  Fx  layer 
echo  trace.  The  unknowns  to  be  solved  for  are  foF  1,  rbFX 
and  rFX. 

[21]  Two  additional  QP  segments  are  now  involved 
where  one  represents  the  joining  layer,  which  smoothly 
joins  the  E  and  Fx  layers  and  the  other  QPS  representing 
the  Fx  layer  from  the  peak  of  the  Fx  layer  down  to  the 
point  where  these  two  layers  are  smoothly  attached.  The 
joining  layer  parameters  depend  on  the  E  and  Fx  layer 
parameters  [see  Dyson  and  Bennett ,  1988].  Thus  three 
QP  segments  are  required  to  determine  the  Fx  layer 
parameters  using  the  inversion  technique  described  here. 
The  equation  for  the  group  path  in  determining  the  Fx 
layer  parameters  contains  the  parts  of  the  raypath  in 
(1)  the  free  space  (region  between  transmitter  and  base  of 
the  ionosphere),  (2)  the  QPS  representing  the  E  layer,  (3) 
the  joining  QPS  layer  joining  the^  and  FI  laws,  and 
(4)  the  QPS  representing  the  Fx  layer. 


[22]  The  equation  for  the  group  path  of  the  propagated  246 
raypathsjbeihg  reflected  from  the  Fx  layer  may  be  written  247 

as  248 


P'  =  -  rl  cos2  (3  -  rQ  sin  (3  +  UE  -  LE 

+  Uj  —  Lj  +  UF 1  —  Lp\  . 


250 


Once  the  Fx  layer  parameters  are  known  the  F2  251 
layer  parameters  may  be  calculated  using  the  same  252 
procedure.  In  all,  up  to  five  QP  segments  are  required  253 
to  determine  the  F2  layer  parameters.  254 


4.  Inversion  of  a  Synthesized  BSI  255 

[24]  To  test  the  inversion  technique  and  its  stability  and  256 
inherent  accuracy,  a  synthesized  backscatter  echo  trace  257 
with  known  QP  layer  parameters  was  used.  The  synthe-  258 
sized  group  path  versus  elevation  angle  echo  traces,  259 


20  40  60 

Elevation  angle  [deg] 


80 


Figure  3.  Simulated  3-D  numerical  ray-tracing  results  using  the  IRI  ionosphere,  where  the  open 
squares  show  the  layer  echo  data  points  chosen  as  inputs  to  the  inversion  technique. 
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Plasma  Frequency  [MHz] 

Figure  4.  Resultant  three-layer  ionospheric  profile  obtained  from  the  inversion  technique  (thick 
curve)  and  the  corresponding  IRI  ionospheric  profile  at  the  transmitter  location  (thin  curve). 


260  shown  in  Figure  1,  were  determined  using  an  analytic 

261  ray-tracing  program  and  the  ionospheric  layer  parame- 

262  ters:  foE  =  3.32  MHz, ymE  =14  km,  and  rE  =  101  +  r°  km 

263  ( E  layer);  foFl  =  4.75  MHz,  ym  F\  -  43  km,  and  = 

264  164  +  ro  km  (Fi  layer);  and  foFl  =  6.45  MHz,  ym  F2  = 

265  68  km,  =  214  +  r°  km  (F2  layer);  where  r°  =  6370  km 

266  (the  radius  of  the  Earth). 

267  [25]  An  operating  frequency  of  12  MHz  was  used  in 

268  the  determination  of  the  backscattered  trace  in  Figure  1. 

269  The  layer  echo  traces  are  clearly  visible  where  the  E 

270  layer  echo  trace  is  the  trace  between  0°  and  12°^leva- 

271  tion,  the  Fx  layer  echo  traces  is  the  trace  between  12°  and 

272  19.7°  and  the  F2  layer  echo  trace ^  between  1J.70  and 

273  29.3°  elevation.  Raypaths  for  this  frequency  that  have 

274  elevation  angles  greater  than  2S^  will  penetrate  the 

275  ionosphere  specified  by  the  above  layer  parameters. 

276  [26]  The  inversion  technique  was  then  applied  to  this 

277  synthesized  data.  Table  1  shows  the  three  data  points 

278  chosen  from  each  of  the  layer  echo  traces  in  Figure  1 

279  which  is  the  minimum  number  of  points  required  for  the 

280  inversion  technique.  The  maximum  elevation  angle  from 

281  each  of  the  layers  is  also  a  required  input.  In  this  initial 


test,  high-accuracy  input  data  were  used  to  validate  the  282 
technique.  283 

[27]  The  results  in  Table  2  clearly  demonstrate  that  the  284 

inversion  technique  converges  correctly  for  the  first  stage  285 
of  the  inversion,  that  is,  the  determination  of  the  E  layer.  286 
After  only  five  iterations  the  technique  has  homed  into  287 
tire  desired  layer  parameters.  The  error  in  group  path  288 
between  the  computed  and  true  values  in  group  path  is  289 
negligible  after  only  five  iterations.  290 

[28]  The  results  in  Table  3  clearly  demonstrate  that  the  291 

inversion  technique  converges  correctly  for  the  second  292 
stage  of  the  processes  determining  the  Fx  layer  param-  293 
eters.  After  only  six  iterations  the  technique  has  homed  294 
into  the  desired  Fx  layer  parameters.  The  difference  in  295 
group  path  after  only  six  iterations  is  negligible.  296 

[29]  The  results  in  Table  4  also  show  that  the  inversion  297 
technique  converges  correctly  in  the  final  stage  of  the  298 
inversion  as  after  only  five  iterations  the  technique  has  299 
homed  into  the  desired  F2  layer  parameters.  The  difference  300 
in  group  path  after  only  five  iterations  is  again  negligible.  301 

[30]  Tables  2-4  clearly  validate  the  inversion  method  by  302 
showing  that  the  new  computed  layer  parameters  converge  303 
closer  to  the  true  values  after  each  iteration.  Figure  2  shows  304 


t5.1 

Table  5.  E  Layer  After  Five  Iterations81 

t5.2 

E,  km 

Error  in  E,  km 

t5.3 

652.8505 

4.81E-7b 

t5.4 

579.0269 

2.81E-7 

t5.5 

568.1541 

3.55E-6 

t5.6 

aFor  the  E  layer,  foE  =  3.55  MHz,  ymE  = 
106.08  km. 

bRead  64.81E-7  as  4.81  x  10“7. 

12.62  km,  and  hmE  = 

Table  6.  F\  Layer  After  Four  Iterationsa  t6.1 


E,  km 

Error  in  E,  km 

t6.2 

966.5592 

-5.59 

t6.3 

953.9551 

-7.85 

t6.4 

955.6044 

3.39 

t6.5 

aFor  the  F\  layer,  foFl  =  5.17  MHz,  ymFl  =  37.74  km,  and  hmFl  = 
196.99  km.  t6.6 
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t7.1  .Table  7.  F2  Layer  After  Nine  Iterationsa 


t7.2 

P',  km 

Error  in  P1,  km 

t7.3 

737.4628 

-1.61E-9 

t7.4 

747.9417 

-2.08E-9 

t7.5 

753.5186 

-1.91E-9 

aFor  the  F2  layer,  foF2  =  8.44  MHz,  ymF2  =  1 12.91  km,  and  hmF2  = 
t7.6  296.56  km. 

bRead  1.61E-9  as  1.61  x  10“9. 

305  the  three-layer  ionospheric  profile  having  the  resultant 

306  layer  parameters. 

307  5.  Inversion  of  Synthesized 

308  Three-Dimensional  Ray-Tracing 

309  Results  Using  the  IRI  Model 

310  [31]  The  backscatter  results  in  Figure  3  show  synthe- 

311  sized  one-hop  sea  scatter  using  a  frequency  of  8  MHz 

312  where  each  point  corresponds  to  a  ray  tube  determined 

313  from  the  three-dimensional  (3-D)  numerical  ray-tracing 

314  program  HIRT  (Homing-In  Ray  Tracing)  [Norman  et  al., 

315  1994].  The  International  Reference  Ionosphere  (IRI) 

316  model  [j Bilitza ,  2001]  was  used  to  provide  a  3-D  repre^ 

317  sentation  of  the  Earth’s  ionosphere.  The  IRI  settin 

318  were:  year  2000;  day  50  (local  summer);  time  1100  LT; 

319  geographic  transmitter  coordinates  (43.4°S,  14j^°E; 

320  55° A)  corresponding  to  the  TIGER  SuperDARN  radar 

321  at  Bruny  Island,  Tasmania,  Australia.  The  synthesized 

322  raypaths  are  determined  in  the  direction  of  the  south 


[32]  Figure  4  compares  the  resultant  three-layer  iono-  325 

spheric  profile  obtained  from  the  inversion  technique  326 
(thicker  line)  with  the  thin  curve  representing  the  IRI  327 
ionospheric  profile  at  the  transmitter  location.  These  328 
results  clearly  show  that  this  new  inversion  technique  329 
produces  accurate  results.  The  results  below  highlight  the  330 
three  sets  of  layer  parameters  representing  the  E ,  F\  and  331 
F2  ionospheric  layers  determined  from  the  inversion  332 
technique  using  the  data  points  represented  by  the  333 
squares  from  Figure  3  and  the  small  number  of  iterations  334 
required  as  well  as  the  high  level  of  numerical  accuracy  335 
obtained  (Tables  5,  6,  and  7)  /  336 

[33]  The  results  in  Figure  5  compare  the  analytic  ray-  337 

tracing  results  obtamed  using  the  three  layer  parameters  338 
determined  from  the  inversion  process  (light  diamond  339 
points)  with  the  3-D  numerical  ray-tracing  results  340 
obtained  using  the  IRI  model  (dark  dashed  points).  341 
The  data  represented  by  the  squares  show  the  values  342 
used  for  the  invasion.  It  is  only  in  the  low-angle  Fx  343 
and  F2  layer  echo  traces  where  there  is  a  significant  344 
difference  in  results.  Ray  tracing  shows  that  there  is  345 
considerably  spreading  of  the  ray  tubes  reflected  at  346 
these  levels  and  this  causing  the  signals  to  be  weak  347 
and  unlikely  to  be  important  for  propagation  experi-  348 
ments  or  communications.  349 


Example  Using  Real  Other  Data 


350 


[34]  The  inversion  technique  was  tested  using  real  351 
backscatter  data  from  the  Tasman  International  Geospace  352 
Environment  Radar  TIGER  [Dyson  and  Devlin ,  2000].  353 


Figure  5.  Resultant  three-layer  analytic  ray-tracing  results  obtained  using  the  parameters 
determined  using  the  inversion  technique  (gray  diamonds)  and  the  corresponding  3-D  numerical 
ray-tracing  results  using  the  IRI  model  (solid  squares).  The  open  squares  represent  the  input  data 
points  chosen  for  the  inversion  process. 
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Figure  6.  TIGER  data  points  for  the  inversion  process  (open  s.q 
analytic  ray-tracing  results  obtained  using  the  parameteg 
(diamonds). 


The  square  data  points  in  Figure  6  represent  the  back- 
scattered  results  received  from  TIGER  and  these  are  the 
data  points  used  in  the  inversion  program  to  produce 
the  following  ionospheric  layer  parameters:  foE 
3.41  MHz,  ym  E  =  17.8  km,  and  rE  =  107.0  + 

(E  layer);  foFl  =6.60  MHz,  ym  F\  =  84.4  km,  and  rFX  = 
197.5  +  ro  km  (FI  layer);  and/oF2  =  6.90  MHz, 

55.9  km,  and  rF2  =  279.5  +  r°  km  (F2 
r°  =  6370  km  (the  radius  of  the  Earth) 

[35]  Figure  7  shows  the  ionospheric  nrofifPtefctained 
by  inversion  and  represented  by  the  above  layer  param¬ 
eters.  The  diamonds  in  Figure  6  represent  analytic  ray¬ 
tracing  results  obtained  using  the  inverted  profile  shown 


L  the  resultant  three-layer 
led  by  the  inversion  technique 

7 


in  Figure  T.  The  data  points  match  up  well  with  the  369 
synthesized  results.  It  is  only  in  the  E  layer  where  there  is  370 
a  difference  between  the  synthesized  layer  echo  trace  and  371 
the  real  data.  372 


Mapping/Modeling  the  Ionosphere 


373 


[36]  The  previous  examples  show  that  a  spherically  374 
stratified  three-layer  model  representation  of  the  iono-  375 
sphere  in  the  region  of  the  radar  can  be  determined  using  376 
the  inversion  technique.  It  will  now  be  shown  how  this  377 
new  inversion  technique  can  be  used  to  determine  a  3-D  378 


4  6 

Plasma  Frequency  [MHz] 


10 


Figure  7.  Resultant  three-layer  ionospheric  profile  obtained  from  the  inversion  of  the  TIGER 
backscattered  data. 


8  of  10 


379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 


xxxxxx 


NORMAN  AND  DYSON:  HF  RADAR  BACKSCATTER  INVERSION  TECHNIQUE 


XXXXXX 


Figure  8.  Raypath  with  an  apogee  height  in  the  Fx 
region  of  the  ionosphere.  The  ray  passes  through  two 
different  E  layers,  represented  by  the  light  gray  and  dark 
gray  grid  point  locations. 


grid  point  model  of  the  ionosphere.  The  inversion 
technique  described  above  assumes  that  the  ionosphere 
is  horizontally  stratified  over  the  region  traversed  by  the 
signals  used  in  the  analysis.  However,  if  horizontal 
gradients  exist,  different  raypaths  will  be  affected  differ¬ 
ently  and  this  will  limit  the  accuracy  of  the  profile 
obtained  by  the  inversion  technique  described  here. 

[37]  On  the  other  hand,  it  is  possible  to  determine  the 
horizontal  ionospheric  gradients  in  the  coverage  area  of 
the  radar  by  applying  the  inversion  technique  to  results 
obtained  at  a  number  of  different  radar  operating  frequen¬ 
cies.  Raypaths  with  different  frequencies  will  be  reflected 
from  different  regions  of  the  ionosphere.  For  example, 
raypaths  with  frequencies  of,  say,  &MHz  will  in  general 
be  reflected  from  the  E  layer  at  ranges  closer  to  the 
transmitter  and  with  higher  elevation  angles  than  raypaths 
having  a  frequency  of,  say,  12  MHz.  By  stepping  through 
the  radar  operating  frequencies  and  determining  the  E 
layer  parameters  at  each  frequency,^  is  possible  using  the 
inversion  technique,  to  determine  the  E  layer  parameters 
at  specified  range  bins,  in  the  direction  of  the  propagated 
signal,  from  the  transmitter  location. 

[38]  Once  having  obtained  a  grid  point  map/model  of 
the  E  layer  in  this  way,  the  Ffla yer  can  then  be  similarly 
mapped/modeled.  The  raypaths  reflected  from  the  Fx 
layer  will  have  passed  through  the  E  region  at  two 
different  ranges,  as  shown  in  Figure  8.  The  Fx  layer 
parameters  can  then  be  determined  using  the  inversion 
technique  and  the  appropriate  two  sets  of  E  region  layer 
parameters  determined  in  the  analysis  of  the  E  region 
data.  By  stepping  through  in  frequency  the  F\  layer 
parameters  can  be  determined  at  each  specified  range 
bin.  Similarly  the  same  process  can  be  applied  for  the 
raypaths  reflected  from  the  F2  layer  where  these  raypaths 
will  have  passed  through  each  of  the  Fx  and  E  regions  at 


two  different  ranges.  The  F2  layer  parameters  can  then  be  414 
determined  using  the  appropriate  two  sets  of  E  region  415 
parameters  and  two  sets  of  Fx  region  parameters.  416 

[39]  Note  that,  as  shown  in  Figures  8  and  9,  the  range  417 

grids  define  specific  regions  in  height  and  range  such  418 
that  many  raypaths  may  travel  through  a  single  grid  419 
component.  The  Segment  Method  Analytical  Ray  Trac-  420 
ing  (SMART)  technique^ Norman  and  Cannon ,  1997,  421 
1999;  Cannon  and  Normqn ,  1997]  which  performs  422 
analytic  ray  tracing  in  a  grid  point  ionosphere  is  ideal  423 
in  this  situation  and  can  be  combined  with  the  inversion  424 
technique  to  produce  a  better  grid  point  ionospheric  map /  425 
model  and  also  to  perform  2-D  analytic  ray  tracing  using  426 
this  new  ionospher^^^pfodel.  427 

[40]  Many  HF  radars  consist  of  multiple  beams  and  428 

this  approach  can  be  used  for  each  beam.  This  means  that  429 
not  only  the  ionospheric  gradients  in  the  direction  of  the  430 
beam  can  be  determined  but  also  the  gradients  at  right  431 
angles  to  the  beam.  The  coverage  area  of  the  radar  can  be  432 
divided  into  a  grid  consisting  of  latitude,  longitude  and  433 
altitude,  where  each  grid  location  will  be  assigned  a  set  434 
of  the  layer  parameters.  Thus  obtaining  a  complete  map  435 
of  the  ionosphere  in  the  coverage  area  is  possible.  The  436 
Segment  Method  Analytical  Ray  Tracing  (SMART)  437 
technique  [Norman  and  Cannon ,  1997,  1999;  Cannon  438 
anchorman,  1997]  which  has  now  been  upgraded  to  439 
perform  3-D  ray  tracing  can  be  used  in  conjunction  with  440 
the  grid  point  ionosphere  to  obtain  realistic  3-D  ray-  441 
tracing  results  and  to  compare  these  results  with  real  442 
radar  data.  443 

[41]  It  should  be  noted  that  several  other  inversion  444 
methods  have  been  proposed  in  the  past;  for  example,  445 
Landeau  et  al.  [1997]  developed  a  technique  known  as  446 
elevation  scan  backscatter  sounding.  Their  approach  was  447 
to  apply  Bayesian  methodology  described  by  Tarantola  448 
[2005]  to  the  inverse  problem,  where  they  applied  it  to  449 
only  a  single  QP  layer  representing  the  F  layer  of  the  450 
ionosphere.  Caratori  and  Goutelard  [1997]  developed  a  451 
relatively  simple  technique  to  determine  the  spatial  iono-  452 
spheric  gradients  however  it  requires  the  addition  of  453 


Figure  9.  Raypath  with  apogee  height  in  the  F2  region. 
The  ray  passes  through  a  number  of  different  E  and  Fx 
grid  point  layer  representations  of  the  ionosphere. 
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vertical  soundings  and  the  inversion  is  limited  to  the  F 
region  of  the  ionosphere. 

8.  Conclusions 

[42]  This  new  inversion  technique  is  robust  in  that  it 
converges  quickly  to  an  accurate  result.  It  has  the  ability 
to  use  many  data  points  from  each  layer  echo  trace  and  in 
general  the  more  data  points  used  the  better  the  inversion 
results  will  be.  The  technique  adds  a  further  dimension  to 
the  analysis  of  HF  radars,  such  as  the  SuperDARN  radar 
network,  where  the  elevation  angle  of  the  received  signal 
is  measured  and  recorded. 
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